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Abstract 

Helicopter loading histories applied to notched metal samples are used as examples, and 
their fatigue lives are calculated by using a simplified version of the local strain approach. This 
simplified method has the advantage that it requires knowing the loading history in only the 
reduced form of ranges and means and number of cycles from the rain-flow cycle counting 
method. The calculated lives compare favorably with test data. 
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Introduction 


The Palmgren-Miner (P-M) theory of damage has long been used to predict the time to crack 
initiation in metals. This rule states that the fatigue failure occurs when the summation of life 
fractions reaches unity. The successful use of this rule requires proper handling of cycle 
counting, overstrain effects, and local notch mean stress effects. In this paper, the P-M rule 
is used, and the above three complexities are included by using the rain-flow cycle counting 
method, by basing the life calculations on data for specimens that have been prestrained [1-3], 
and by using the local strain approach [4-6], respectively. The local strain approach focuses 
attention on the stresses and strains that occur locally at a stress raiser of interest, and the 
S-N curve used is a strain versus life curve. 

The rain-flow cycle counting method is a procedure for Interpreting an irregular load versus 
time history as a collection of events {called cycles) to which fatigue damage can b^ assigned. 
In this method [7], cycles are counted depending on the comparison of two adjacent ranges 
as illustrated in Fig. 1, which also defines the range and mean of a cycle. If the first range is 
less than or equal to the second, a cycle is counted and the corresponding peak and valley 
are discarded for the purposes of further cycle counting. 

Figure 2 illustrates this process for a simple loading history. First, the history is reordered to 
start with the highest peak or lowest valley as in (b). Cycle counting then proceeds by moving 
forward in the history. If a cycle is counted, as in (c), its range and mean are recorded and its 
peak and valley are removed from the history. This procedure, as described in detail in Ref. 
[7], always yields a major cycle between the highest peak and lowest valley, and smaller cy- 
cles temporarily interrupt larger ones. 

By using the rain-flow cycle counting method, a service loading history can be reduced to a 
convenient compact form. The compact description is in the form of a matrix giving range, and 
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mean, and numbers of rain-flow cycle. Table 1 illustrates such a matrix for the loading history 
of Fig. 2. As a second example, a portion (about one-third) of an actual helicopter loading 
history is shown in Fig. 3, and the resulting rain-flow matrix in Fig. 4. 

As noted above, an area of importance in predicting fatigue life is the overstrain effect. This 
effect, which is caused by the higher stress levels, needs to be considered since it increases 
the damage done by the lower stress levels [1,2]. Figure 5 illustrates this effect for titanium 
6AI-4V. The lower curve represents test data for specimens which have been plastically 
strained, but only to a life fraction of a few percent. There is nevertheless a large effect on the 
life for subsequent testing at a lower level due to the prestrain causing damage to the material 
at the microstructural level. 

A known load versus time history is necessary for analysis of fatigue life using the local strain 
approach. But the history may be lengthy, and there are no restrictions on the degree of ir- 
regularity of the time variation, The local strain approach predicts crack initiation life and 
assumes that fatigue life is controlled primarily by the local notch surface strain and mean 
stress, not the nominal (average) stress. Emphasis on local notch behavior Is crucial as this 
permits rational analysis of local notch yielding and its effect on the local notch mean stress. 
This method in its complete form requires knowing the loading history in full length. 

A simplified method for calculating fatigue crack initiation life based on the local strain ap- 
proach can be used [5,8]. This method has two very distinct advantages. First, only a rain-flow 
matrix in the compact form of range-mean values is required as the input information. Al- 
though some detail is lost, such a matrix can be used with the local strain approach to place 
upper and lower bounds on the life that would result from the analysis of the original, un- 
summarized history. This principle is explained in detail in Refs. [3,9]. Note that this compact 
form is much easier to handle and store than the full history in the form of a time sequence 
of peaks and valleys. Secondly, the life calculations are simpler and more economical. 
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BC< AB BC > AB 

No cycle AB = cycle 


For cycle A - B 

Peak = S A 

Valley s S B 

Range = AS =S A - S B 

Mean = So = (S A +S B )/2 


Figure 1. Condition for recording an event during rain-flow cycle counting 
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Load Units Load Units Load Units 



Figure 2. Example of rain-flow cycle counting from the ASTM standards |7]: Before the cycle 

counting begins, the most extreme point in the history should be located, and the history 
arranged to start and Finish at this point as shown in (b). 
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Figure 3. Portion of the modified maneuver history 
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Table 1. Rain-flow matrix for the loading history of Fig. 2 (7] 
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N , Cycles 


Figure 5. Strain vs. life test data and curves for TI-6AI-4V 
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In this paper, two helicopters loading histories are chosen for analysis, and their fatigue lives 
are calculated using the simplified method. Also, the calculated lives are compared with test 
data. All life calculations are based on a strain-life curve, which represents failure of small 
unnotched axial test specimens. Lives so predicted for notched members correspond to initi- 
ation of an easily detectable crack. In general, the remaining life for crack growth also needs 
to be calculated, but this was relatively short and so was neglected for the particular speci- 
mens analyzed here. 


Life Calculation Methodology 


The simplified version of the local strain approach used assumes that fatigue life is controlled 
primarily by the notch surface strain, and it considers plasticity and mean stress effects in a 
fairly complete manner. Note that the traditional S-N approach crudely handles plasticity, and 
also neglects the special mean stress effects which arise from loading sequence. Life calcu- 
lations by the local strain approach consists of two steps. First, the local notch stress and 
strain histories must be predicted, and second the life corresponding to the local stress and 
strain histories must be estimated. 


Local Strain Approach 

Figure 6 illustrates the initial and most difficult step of estimating the strain-stress history. 
This step is difficult because it requires specific handling of the complex nonlinearity relating 
load, strain, and stress. In order to achieve the above task, a cyclic stress curve [10] for the 
material is needed: 
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( 1 ) 


e* = 



where a, are amplitudes of strain and stress, respectively, E is the elastic modulus, and A 
and s are material constants. Next, by employing Eq.1 and with the aid of NeubeKs rule, a 
curve relating nominal stress, S, and the local notch stress and strain Is obtained: 


(*<s a ) 2 

ff a c a= £ 


( 2 ) 


where k, is the elastic stress concentration factor. 

Figure 6(b) shows these two curves. These curves are then used to estimate the local 
stress-strain response at the notch by following the loading history while modeling the 
hysteresis looping behavior of the material. For the example of Fig. 6, the irregular loading 
history of (c) results in S versus e and a versus e as shown in (d) and (e). Note that there is 
a set of closed hysteresis loops, such as 2-3-2', 6-7-6', 5-8-5', and 1-4-1' for this example. Each 
such loop is identified as a cycle, and the cycles so defined are the same as would be ob- 
tained from applying rain-flow cycle counting to the load (S) versus time history. 


Each cycle now has a known strain range, As = 2e„ and mean stress, a 0 , as shown for cycle 
6-7-6' in Fig. 6 (e). The life, N, corresponding to each combination of e, ,o 0 can be obtained from 
a strain-life curve [ 10]: 


E a = (2N') b + e/(2A/') c (3) 

where e, is the strain amplitude corresponding to a closed loop, N' is the life in cycles for zero 
mean stress, and a,', b, t{ and c are additional constants for the material. If the rule of Morrow 
(11] is used to account for mean stress effects, the fatigue life as adjusted for mean stress <r 0 
can be estimated by. 
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(4) 


N = N'( 1 



where N is the final adjusted life. 

The final step is then to apply the P-M rule. For a loading history that is assumed to repeat 
until failure occurs, the rule takes the form: 



L per block J 


(5) 


where n, is the number of occurrences per block of a cycle corresponding to life N„ and B is 
the unknown number of blocks (repetitions) to failure for the irregular history. 


Simplified Method 


As mentioned earlier, the local strain approach generally requires a knowledge of the full 
sequence of the loading history. However, in the simplified method, the loading history may 
be used in concise matrix form. The first step is then to determine the rain-flow matrix that 
contains the information on range and mean loads of rain-flow cycles. Using this information, 
upper and lower bounds on life can be calculated [3,9). As long as these bounds are reason- 
ably tight, which will be the case in most practical situations, the more detailed simulations 
as in Fig. 6 are unnecessary. 

The principle behind this bounding is shown in Fig. 7 for one cycle, namely 6-7 from Fig. 6. The 
guiding principle is that both load-strain loop 6-7 and stress-strain loop 6-7 must lie within the 
corresponding loop for the largest cycle in the history, namely 1-4. The load limits S 6 and S 7 
are known; therefore limits can be placed on the mean strain of cycle 6-7. As shown in Fig. 
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Figure 6. Illustration of local-strain approach for an Irregular history: Notched member (a), 

' having cyclic stress-strain and load-strain curves as in (b), is subjected to load history 

(c). The resulting load-strain response is shown in (d). and (e) is the local- notch stress 
strain response [6], 
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7(a), loop 6-7 can be so far to the right that it is attached at A, or so far to the left that is at- 
tached at B. Similarly, the same line of reasoning can be applied to a load-stress loop as in 
Fig. 7(c), where loop 6-7 could be so low that it is attached at A, or so high that is attached 
at B. Figure 7(b) shows the extreme stress-strain loops which satisfy both sets of constraints, 
so that these correspond to the upper and lower bounds on local notch mean stresses for 
cycle 6-7, and o M . 

Knowing the upper and lower bounds on the mean stress for cycle 6-7 of the example allows 
bounds to be placed on the life, N, from Eq. 4. The upper bound on N is similarly obtained for 
all cycles in the history, and these are used with the P-M rule in the following form of Eq. 5 to 
obtain the upper bound on life for the irregular history. The same procedure, but using the 
lower bound N's for each cycle, can be used to obtain the lower bound on life for the irregular 
history. The above procedure is explained in detail in Refs. [3,9], 

If the local notch plastic strains are small during the history, the stress-strain loop for the 
major cycle, such as 1-4 in the example, is reduced to a straight line, and therefore the upper 
and lower bounds on life are then identical. Also, at high load levels, the cycles causing most 
of the damage may not have significant mean stresses due to the large plastic strains, and the 
bounds will then be close. Hence, the widest separation between the bounds is expected at 
intermediate load levels. Note that if most of the damage is done by low-level cycles, then the 
degree of separation of the bounds will be greatest. Also, if all but a negligible fraction of the 
damage is done by the major cycle , then the two bounds are again identical [3,9]. 
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Figure 7. Illustration of placing bounds on the mean stress of a subcycle: This illustration is 

based on Fig, 6, Note that the sequence of the applied loads is not known, and the mean 
stress for 6-7 must lie between the a 0A and <7,^ values shown. 
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Material, Specimens and Loading Histories Analyzed 


Two helicopter loading histories are analyzed. Fatigue life data for both Helix and a severe 
maneuver history, applied to plate-with-hole specimens having elastic stress concentration 
factors between 2.4 and 3.92, are available in Ref. [12] and from tests done at the University 
of Dayton Research Institute (UDRI). These data include several levels of maximum nominal 
stress, for titanium 6AI-4V. 

Two loading histories, Helix and maneuver, are used in fatigue life calculations. Helix is a 
standard helicopter loading spectrum obtained from Ref. [12], and the maneuver history was 
obtained by the University of Dayton Research Institute (UDRI). These two histories are further 
explained below. 


Materials Properties Used 


Figures 5 and 8 show strain vs. life and cyclic stress vs. strain data for the titanium 6AI-4V 
material used, and also curves corresponding to Eqs. 1 or 3 fitted to the data. Constants cor- 
responding to these curves are given in Table 2 for constant amplitude tests and also for tests 
on prestrained material. Some constant amplitude tests were done at Virginia Tech on the 
same material. These data are also shown in Figs. 5 and 8, and they agree well with the data 
obtained from Ref. [13] and also with the fitted constants. In the case of the prestrained ma- 
terial, the data and constants were obtained from Ref. [14]. 
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Figure 8. Cyclic stress-strain test data and curves for TI-6AI-4V 


Fatigue Life Estimates for Helicopter Loading Spectra 


16 


, Stress Amplitude , ksi 




Table 2. Cyclic stress-strain and strain-life constants for Ti 6AI-4V 


Symbols, Units 

Titanium 6AI-4V 


Const. Ampl. 

Prestrained 

E, GPa (ksi) 

113.8 (16500) 

113.8(16500) 

A, MPa (ksi) 

1327 (192.4) 

1702 (246.9) 

s 

0.0755 

0.127 

e', 

6.22 

2.802 

c 

-1.01 

-0.860 

a'„ MPa (ksi) 

1523 (220.9) 

2207 (320.0) 

b 

-0.0763 

-0.126 
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General Description of Helix 


Helix [12] is a standard loading sequence for the main rotors of helicopters with articulated 
rotors. Helix represents a loading history for a 190.5-hour (2,132,024 cycle) sequence of 140 
flights. Each flight in the sequence represents one of either training, transport, antisubmarine 
warfare, or search and rescue. Each of these appears in the sequence in three different 
lengths of 0.75 hour, 2.25 hours, and 3.75 hours. There are twelve unique flights, which are 
applied in a specific number of repetitions and sequence to obtain the total 140 flights. Figure 
9 shows the load vs. time history for portions of a transport flight in Helix. 

Helix is composed of 24 unique maneuvers, which are repeated in various sequences and 
numbers of repetitions to compose the various flights. The maneuvers such as take off, for- 
ward flight of various load levels, turns, etc. each consists of a mean level and a relatively 
small number of cycles. These cycles occur at one or more stress amplitudes, with the num- 
ber of cycles being between 1 and 40. 

Helix has a relatively high mean levels for the various maneuvers, these mostly ranging from 
60 % to 68 % of the maximum nominal stress in the spectrum , . Helix reaches 100 % of 

the S^,, level at least once in each flight and returns to -20 % at the end of each flight. Hence, 
Helix has a large ground-air-ground cycle and a large number of cycles at relatively high 
mean levels. Table 3 gives the range-mean matrix for Helix from rain-flow cycle counting. The 
matrix entries in Table 3 were obtained from all of Helix by dividing each by 140, the number 
of flights. These range and mean values correspond to the history scaled so that S m<> = 100 
units. 
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Time 


Figure 9. Example of the loading history for portions of a transport flight in Helix [12] 
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Table 3. Range-mean matrix for Helix from rain-flow cycle counting 
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General Description of Maneuver History 


A loading history for the tail rotor pitch beam of an AUH-76 helicopter was selected as repre- 
sentative, and loading histories from each of 30 distinct severe maneuvers were assumed to 
occur once each in a specific sequence. This produced a loading history containing 33,470 
cycles, which was then modified by UDRI to eliminate minor events, shortening it to 8777 cy- 
cles, that is, 8777 peaks and 8777 valleys. This new history Is called the modified maneuver 
history. Figure 3 shows about one third of the modified maneuver history. Figure 4 shows the 
complete range and mean matrix from the rain-flow cycle counting. Note that the maneuver 
history is scaled so that the highest load is 1 unit, which results In the lowest load being -0.516 
units; therefore, the overall range is 1.516 units. 

Fatigue life calculations were done for both the original and modified histories for titanium 
6AI-4V. It was found that the differences are small, and the two histories appear to cause the 
same damage; therefore, all further testing and analysis is based on the modified history. 


Upper-Lower Bound Analysis of Helix 


Fatigue lives for the plate-with-hole specimens of titanium 6AI-4V subjected to Helix were 
calculated using the simplified {upper-lower bound) version of the local strain approach. As 
already discussed, only the rain-now matrix of the history is required as input information for 

fatigue life analysis. 

Rather than analyzing all of Helix as a single loading history, it is expedient to analyze each 
of the twelve unique flights separately. Therefore, the first step is to determine the rain-How 
matrix of each unique night. Then, by using these matrices . upper and lower bounds can be 
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determined for each unique flight. These lives then combine to obtain overall bounds by con- 
sidering the number of times each unique flight is repeated in Helix. Note that this combining 
of flights works because all of the flights return to the -20% level. This gives a common mini- 
mum stresses and strains locally at the notch for each flight, which results in no sequence 
effects from one flight affecting another. 

Another option is also available for predicting the lives. In this option only one matrix is used 
which is derived from all of Helix. Comparison between these two options shows no signif- 
icant difference; therefore, the simpler one of a single matrix was adopted. Hence, the rain- 
flow matrix for an average flight in the form of Table 3 is used for the purpose of calculating 
the upper and lower bounds on life. 

The resulting calculated lives are plotted in Fig. 10 and are given in Table 4. Note that one 
curve is based on constant amplitude strain-life data, and the other on prestrained data. Fig- 
ure 10 plots k, S m „ so that the data for various k t can be shown on the same plot. This is ex- 
pected to be valid based on Neuber's rule, Eq. 2, as long as net section yielding does not 
occur, in which case Eq. 2 is not valid. 

The data obtained from RAE scatter over a broad range, and the lives tend to be longer than 
those obtained from UDRI. Also, from Fig. 10, general agreement is obtained between data 

and analysis. 


Upper and Lower Bound Analysis of Maneuver History 


Figure 11 shows the calculated results as well as test data. Table 5 gives the calculated lives 
for notched specimens (k,= 2.5). Reasonable agreement is obtained by comparing the test 
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MPa 



Figure 10. Analysis of Helix compared to test data for Ti 6AMV: The lines shown are middles 
of bounds. 
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Table 4. Calculated flights to failure for Helix for notched specimens (k t =2.5) 


S Mr MPa (ksi) 
(net area) 

Helix 

Lower Bd. Upper Bd. 

(a)Ti 6AI-4V, Const. Ampl. Strain-Life Curve 

317 (46) 

26193 

26213 

358 (52) 

4415 

4433 

407 (59) 

857 

876 

455 (66) 

221 

240 

517 (75) 

46.9 

62.8 

(b) Ti 6AI-4V, Prestrained Strain-Life Curve 

317 (46) 

681 

687 

455 (66) 

32.7 

36.0 

517 (75) 

11.8 

13.9 

530 (77) 

9.59 

11.5 


Fatigue Life Estimates for Helicopter Loading Spectra 


24 




data and the calculated values, especially when prestrained data are used. An exception is 
at the lowest stresses where these calculations are conservative. 

Figure 12 shows the distribution of numbers of cycles vs. ranges. During the determination 
of the fatigue lives, it was noted that most of the fatigue damage was done by the higher range 
levels, and none by the lower levels. In this regard, it was also noted that a potentially large 
saving of test time can be realized by eliminating lower level non-damaging stress cycles from 
the load history. Therefore, a rain-flow filtering was done on the maneuver history. 

The basis of this filtering was the fatigue "damage", more properly called the usage fraction, 
which is obtained from the P-M rule calculation. Usage fraction is defined as: 

N, 

usage fraction — — (6) 



where n, is the number of cycle applied at a stress level corresponding to life N„ and I is the 
number of different discrete stress levels. Figure 13 shows the usage fraction at each level vs. 
range for the modified history. 

As is evident from Figs 12 and 13, most cycles occur at the lower ranges, whereas most 
damage is estimated to be done at higher ranges. Therefore the modified history was short- 
ened by filtering all the ranges less than 0.45 units, that is, 30 % of the largest rain-flow range 
of 1.516 units. This filtered history contains 510 cycles, that is, 510 peaks and 510 valleys. 
Figure 14 shows the filtered history. By filtering the lowest stress levels of the modified history, 
the usage fraction is decreased only slightly, specifically by about 1 %. 

Note that both modified and filtered histories produce similar fatigue lives, since most of the 
damage is done at higher levels that are not filtered out in the latter. Also, from Fig. 11 com- 


Fatigue Life Estimates for Helicopter Loading Spectra 


25 



MPa 



Figure 11. 


Analysis of maneuver history compared to test data for H 6AI-4V: The lines shown 
are middles of bounds. 


Fatigue Life Estimates for Helicopter Loading Spectra 


26 






RANGE 





29 


parison of the test data for modified and filtered history at one stress level, where both were 
tested, shows quite good agreement, which indicates that the filtering was a success. 


Discussion 


The Helix spectra is typical of a case where the bounds would tend to be relatively widely 
separated. This is due to most of damage being done at a relatively low level within the 
spectrum. However, even in this unfavorable situation, the bounds are reasonably tight for 
both Helix and maneuver history. For real irregular loading histories , wide separation of 
bounds is not likely to very often be a difficulty with this approach. 

From Figs. 10 and 11, it is seen that the calculated lives and the test data are In good agree- 
ment when the prestrained data are used. The only exceptions are in the lowest levels, where 
the data approach or even exceed the calculation based on the constant amplitude strain-life 
data. The reason for this trend appears to be that these stress levels are so low that even the 
highest stress level (ground-air-ground) is approaching the endurance limit, and no prestrain 
effect due to this highest level occurs. The interpretation is made that when the major cycle 
is sufficiently low, an endurance limit is expected for the spectrum loading. Based on the 
above discussion, small cycles with amplitudes below the endurance limit can cause fatigue 
damage if preceded by a major cycle substantially above the endurance limit. 

The calculated lives tend to be reasonably accurate at high stress levels, although there is 
considerable scatter. For Helix on titanium, the highest stress levels used were limited due to 
static failure, even though the most damaging level was still relatively low. This contrasts with 
the situation for the maneuver history, where the most damaging levels were calculated to 
be those near the maximum rain-flow range in the history, (see Fig. 13) 
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Table 5. Calculated repetitiona to failure for maneuver hietory for notched specimens (k t =2.5) 


S m „, MPa (ksi) 
(net area) 

Modified Sequence 
Lower Bd. Upper Bd. 

(a)Ti 6AI-4V, Const. Ampl. Strain-Life Curve 

276 (40) 

67700 

67700 

327 (47) 

7940 

7940 

379 (55) 

1410 

1410 

455 (66) 

291 

315 

586 (85) 

59.3 

68.9 

827 (120) 

6.14 

13.8 

(b) Ti 6AI-4V, Prestrained Strain-Life Curve 

276 (40) 

8770 

8930 

379 (55) 

888 

934 

827 (127) 

7.89 

14.1 
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Similar comparisons from previous work [3] for Helix applied to notch specimens of an alu- 
minum alloy are shown in Fig. 15. The trends are the same as for titanium, reinforcing the 
comments just made. 

The above discussion and comparisons between test data and calculated lives suggest that 
the upper and lower bound approach represents a useful method with distinct advantages for 
predicting fatigue crack initiation. 


Conclusions and Recommendations 


The following conclusions and recommendations are drawn based on the above analysis and 
discussion: 

1. Overstrain effects caused by the higher stress levels in a load spectrum need to be con- 
sidered since these increase damage at the lower stress levels. This is especially true for 
cycles below the endurance limit where some cycles in the load spectrum are above the 
endurance limit. 

2. The simplified version of the local strain approach (upper and lower bounds) should be 
used more widely, since only a rain-flow matrix is required as input information, and since 
it is easy to program on a digital computer. This approach is more economical than 
analysis of the large amount of data involved in full time sequence load histories. Based 
on the comparison of test data and calculated lives, the accuracy of the simplified method 
was either reasonably accurate or conservative depending on the level. 
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Smoi i Maximum Nominal Strass , MPa 



Figure 15. Analysis [3] of Helix compared to test data [12] for 2024-T4 aluminum notched speci- 
mens with k =2.5: The curves represent full analysis based on local strains, and the 

bounds from the simplified method are also shown. 
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ABSTRACT 


A computer program (UPLO) for predicting upper/lower bounds on life 
is provided. The program takes the result of rainflow counting in the 
form of a matrix, that is, given numbers of cycles at each combination 
of range and mean values, and then uses this information by appling a 
local strain approach method to place upper and lower bounds on life. 

The input values are defined, and one example problem is attached. 
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INTRODUCTION 


Consider a notch member as in Figure A- la, subjected to an 
irregular variation of nominal stress, S, with time. The goal is to 
predict the upper and lower bounds on life. A simplified version of 
local strain approach described In Reference A-l was used to predict the 
upper and lower bounds on life. The first step is to summarize a 
lengthy history, using the rain-flow counting method, into a compact 
form of a matrix giving combinations of range and mean values. This 
step must be done using a rain-flow computer (RAINF) program which feeds 
its output to the upper and lower bounds program as data. The above 
matrix is then used with the local strain approach to place upper and 
lower bounds on life that would result from the analysis of the original 
unsummarized history. The principle behind this bounding is illustrated 
by Figure A-2 with the aid of Figure A-l. 

The computer program explained in this manual is the same program 
that was previously used in Reference A-2. 
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Figure A-l - Illustration of local strain approach for an Irregular 
load vs. time history. Notched member (a), having cyclic stress-strain 
and load-strain curves as In (b). Is subjected to load history (c). The 
resulting load-strain response Is shown In (d) and the local notch 
stress-strain response In (e). 
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Figure A-2 - Illustration based on Fig. 1 of placing bounds on the 
mean stress of a subcycle when the sequence of the applied loads Is not 
known. The mean stress for cycle 6-7 must lie between the values of 
and 


'oA 


oB 
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DETAILS OF EQUATIONS 

The equations associated with the program are listed below to 


provide necessary information for the user. 


1. Cyclic stress-strain curve: 


_ °a . A* 1/s 
£ a - r + ( A > 


( 1 ) 


where 

a = stress amplitude 
d 

E = strain amplitude 
a 

E = elastic modulus 
A = cyclic strength coefficient 
s = cyclic strain hardening exponent 

2. General strain-life curve: 

e, . ijf (2N*) b ♦ e ' f (2N*) c ( 2 ) 

where 

o' = fatigue strength coefficient 

b = fatigue strength exponent 

E ' f = fatigue ductility coefficient 

c = fatigue ductility exponent 

N* = life in cycles for zero mean stress 

3. Life equation due to Morrow, considering the effect of mean stress. 


N = N* (1 - -t - ) 1/b 

0 f 


( 3 ) 


where 

N = final life 

0 = mean stress 

o 
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4. Palmgren-Miner rule: 


B l TT “ 1 (4) 

per block N 1 

where 

= number of occurrences of a cycle corresponding to life 
B = unknown number of blocks to failure 


5. Neuber's rule: ? 

< k t V 

Va " E 

where 


( 5 ) 


k^ = stress concentration factor 
S a = nominal stress amplitude 
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PROGRAM PROCEDURE OF UPPER/LOWER BOUND CALCULATIONS 


The detailed procedure of the upper/lower program is given below. 
The example of cycle 6-7-6' of Figure A-l is further employed as an 
example with the aid of Figure A-2. 

It is convenient to write Equations 1 and 5 in general form without 

subscripts: 


e = o/E = (a/A) 1/S 


( 6 ) 


ae 


(k t S ) 2 

E 


( 7 ) 


The values of the constants E, A, s and k^ are of course unchanged. 

Combining Equations 6 and 7 gives a relationship involving only 
strain, e, and nominal stress, S. 


Z 




(k s ) 
EeA 


2 1/s 


(8) 


Considering Figure A-2, the goal is to determine the bounds on mean 
stress, such as a QA and o 0 g. Point 1 corresponds to the maximum load in 
the history and 4 to the minimum load in the history. As a convenience, 
it is assumed that the largest absolute value of load is positive. If 
not, then what follows will need to be modified with appropriate sign 
changes. Note that the load history is known, which implies that S 
values are known for all calculations, so that the unknowns are the 
a and e values. These calculations take advantage of the fact that 
various loop curves in either Figure A-2a or b have the same shape, 
which is that of the corresponding curve from Figure A-l expanded with a 
scale factor of two. 
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To obtain the unknowns for point 1 let 


s = s l 

e = e l 

a = 


(9a) 

(9b) 

(9c) 

Substitute Equation 9 into Equation 6 and 8, and solve for from 

Equation 8. Then using Equation 6, solve for oj . Equations 6 and 8 
are solved by Newton's method, as direct solutions are not possible. To 
analyze the range of major cycle 1-4-1', let 

<. aS 1-4 

* ' 2 


(10a) 

Ae l-4 
e " 2 


(10b) 

Ao 1-4 
o - 2 


(10c) 

Using a parallel procedure to that just described, Ae^ and 
obtained. Then the stress and strain at point 4 are 

Ao 1_4 are 

c 4 " e l " Ae l-4 


(11a) 

°4 = °1 " Ao 1-4 


(11b) 

To analyze the range of minor cycles, such as 

6-7, let 


_ aS 6-7 

^ - O 


(12a) 
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Ae 6-7 

e - 2 

(12b) 

* 

Ao 6-7 
0 - 2 

(12c) 


Using the same procedure, Aeg _-j and Aog _-j are determined. 

Once the stress and strain at points 1 and 4 are obtained, the 
points of attachment of loops A and B in Figure A -2 must be determined. 
In order to determine point of attachment of loop A, let 


, S 1- S A 
b - 2 

(13a) 


e l ' e A 

e - 2 

(13b) 


CT 1 - a A 

a - 2 

(13c) 


Then substitute Equation 13 into Equations 
o A . Then to find the point of attachment 

6 and 8 and obtain and 

of loop B, let 


r S B ' S 4 
b - 2 

(14a) 


E B " e 4 
£ - 2 

(14b) 

-- 

°B " a 4 
o - 2 

(14c) 


Again substituting into Equations 6 and 8, 

obtain eg and og . 
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The bounds on mean stresses are then 


a °6-7 

a oB " °B ' 2 


(15a) 


°oA ' °A 



(15b) 


Next, into Eq. 2 substitute 


t 


a 



(16) 


and obtain N, the life for zero mean stress. Finally, substitute this 
N* and o Q g into Equation 3 to obtain the lower bound in life, N, for 
cycle 6-7-6'. Similarly, substitute N* and a QA to get the upper bound 
on N. 

Following a similar procedure for all cycles smaller than the major 
one then allows the P-M rule. Equation 4, to be employed once with all 
of the lower bound N values, and a second time with all of the upper 
bound N values, to obtain bounds on the calculated number of blocks 
(repetitions) to failure, B. 
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DEFINITION OF INPUT DATA 


Line No. 

Of Read Variable 

Statement Name Explanation Comment 


3 

4 

5 


6 


7 


9 


L (Number of columns)x( number of rows) 

in the range/mean matrix 

KT Stress concentration factor, kj. 

AMAX Load scale factor 


RMMIN Lowest mean value in the matrix 

LEVLM Number of columns (mean values) 

in the range/mean matrix 

XIM Constant increment between mean values 

RMIN Lowest range value in the matrix 

LEVLR Number of rows (range values) in the 

range/mean matrix 

XIR Constant increment between range values 
EM E 


Factor which mul- 
tiplies by the 
range and mean 
values from the 
matrix to give 
nominal stresses, 
S in ksl , defined 
consistently with 


All these values 
should be the same 
as the values in 
the range/mean 
rainflow matrix 
used as input. 


ksi 


A A 


ksi 


SH s 

SFP a ' f ksi 

BS b 


EFP e' f 

CS C 

MM ( ) Elements of the matrix in one 
dimensional form 


For Example: If a 

matrix has 10 
columns, then 
MM(1) = A( I , 1) 
MM(2) = A(l,2) 
MM(10) = A( 1 , 10) 
MM ( 1 1 ) = A(2, 1) , 
etc. 

o^, and AMAX. All 


Note: Other consistent stress units may replace ksi for E, A, 

other inputs are dimensionless. 
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EXAMPLE PROBLEM 


A rainflow matrix of combinations of range and mean values for the 
standard helicopter load spectrum Helix [ A-3 ] is used to determine the 
upper and lower bounds on life. Table A-l shows this matrix which 
contains 30 rows and 9 columns. Table A-2 shows the input for this 
example. The entire program listing and output are also attached as 
Table A-3. The output lists the upper bound mean stresses, such as 
a oB ™ ^9 ure A_ 2b and also the lower bound mean stresses, such as 
°oA figure A- 2b. It also gives the percentage of the total "damage," 
that is, the Palmgren-Miner usage factor, corresponding to the upper and 
lower bound mean stress cases for each particular combination of nominal 
stress range and mean. 

The load scale factor AMAX is selected to scale the load history to 
any desired magnitude. To select a value, it is necessary to know how 
the original rainflow matrix used as input Is scaled. AMAX is 
specifically the constant which the values in the input matrix are 
multiplied by to obtain the desired nominal stresses, S. For this 
example, the largest cycle has a mean value from the matrix of 0.40 and 
a range of 1.20. Hence, its maximum is 1.00 and its minimum is -0.20. 

If a load scale factor of AMAX = 59 is chosen, this causes the analysis 
to be done for S max = 59x1.00 = 59.0 ksi, which corresponds to the 
largest range in the load history being a S = 59x1.20 * 70.8 ksi. 
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Table A-l. Range-Mean Matrix for Helix from Rain-Flow Cycle Counting 
(Cycles per Flight, Average)* 


Mean 


Range 

.04 

.08 

.12 

.16 

.20 

.24 

.28 

.32 

.36 

.40 

.44 

.48 

.52 

.56 

.60 

.64 

.68 

.72 

.76 

.80 

.84 

.88 

.92 

.96 

.100 

1.04 

1.08 

1.12 

1.16 

1.20 

ALL 


.40 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


■ 44 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


,48 

.52 

.56 

.60 

.64 

.68 

.72 

ALL 

0 

1 

0 

2 

16 

2 

0 

21 

0 

0 

1 

0 

4 

0 

0 

5 

0 

0 

0 

0 

2 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

2 

0 

43 

30 

1742 

1348 

27 

223 

3425 

0 

1 

0 

1 

0 

1 

0 

3 

0 

14 

5 

453 

2613 

155 

12 

3252 

o 

0 

0 

5 

20 

1 

0 

26 

0 

24 

6 

65 

7785 

460 

4 

8344 

0 

0 

1 

0 

15 

0 

0 

16 

0 

6 

8 

30 

6 

24 

0 

74 

0 

1 

1 

0 

0 

0 

0 

2 

0 

10 

4 

28 

0 

0 

0 

42 

0 

1 

1 

1 

0 

0 

0 

3 

0 

1 

1 

5 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


12 


103 


59 


2333 11811 


671 


239 


15229 


*The matrix entries were obtained from those for all of Helix by dividing each by 140, the 
number of flights. 
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Table A -2. Input for Example Problem 


270 

2.5 

59 . 

0 . 4 , 9 * 0 . 04 , 0 . 04 , 30 * 0.04 

15838 .. 192 . 4 . . 0755 . 220 . 9 , -. 0763 , 6 . 216 , - 1.0101 
0 , 0 , 0 , 1 , 0 , 2 , 16 , 2 , 0 , 

0 , 0 , 0 , 0 , 1 , 0 , 4 , 2 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 2 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 2 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 1 , 0 , 

0 , 12 , 0 , 43 , 30 , 1742 , 1348 , 27 , 223 , 

0 , 0 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 

0 , 0 , 0 , 14 , 5 , 453 , 2613 , 155 , 12 , 

0 , 0 , 0 , 0 , 0 , 5 , 20 , 1 , 0 , 

0 , 0 , 0 , 24 , 6 , 65 , 7785 , 460 , 4 , 

0 , 0 , 0 , 0 , 1 , 0 , 15 , 0 , 0 , 

0 , 0 , 0 , 6 , 8 , 30 , 6 , 24 , 0 , 

0 , 0 , 0 , 1 , 1 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 10 , 4 , 28 , 0 , 0 , 0 , 

0 , 0 , 0 , 1 , 1 , 1 , 0 , 0 , 0 , 

0 , 0 , 0 , 1 , 1 , 5 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

1 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 , 
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Table A-3. Program Listing and Output for Example 


l 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


CtJOB 

C 

C 

C 


A1336 ,P*100 

UPPER/LOWER LIFE PREDICTION PROGRAM (UPLO) 


c 

INPUT 



c 

c 

DATA 

LINE 

1 . 

L=NO OF COLUMNSXNO OF ROWS IN THE RANGE/ 
MEAN MATRIX. 

c 

DATA 

LINE 

2 . 

KT=STRES5 CONCENTRATION FACTOR. 

c 

DATA 

LINE 

3. 

AMAX=LOAD SCALE FACTOR .(KSI) 

c 

c 

c 

c 

DATA 

LINE 

4. 

RMMIN=LOWEST MEAN VALUE IN THE MATRIX. 
LEVLM=NO OF COLUMNSCMEAN VALUES) IN THE 
RANGE/MEAN MATRIX. 

XIM*CONSTANT INCREMENT BETWEEN MEAN VALUES 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


THE 

IN 

DATA LINE 


DATA LINE 5. 


RMIN-LOMEST RANGE VALUE IN THE MATRIX. 

LEVLR=N0 OF ROWS(RANGE VALUESJIN THE RANGE/ 
MEAN MATRIX. 

XIR=CONSTANT INCREMENT BETWEEN RANGE VALUES. 
ABOVE VALUES SHOULD BE THE SAME AS THE VALUES 
THE RANGE/MEAN RAINFLOW MATRIX USED AS INPUT. 

4. EM=ELASTIC MODULUS(KSI) 

A=CYCLIC STRENGTH COEFF(KSI) 

SH s CYCLIC STRAIN HARDENING EXPONENT 
SFP=FATIGUE STRENGTH COEFF(KSI) 

BS=FATIGUE STRENGTH EXPONENT 
EFP=FATIGUE DUCTILITY COEFF 
CS=FATIGUE DUCTILITY EXPONENT 
MM( )=ELEMENTS OF MATRIX IN ONE DIMEN- 
SIONAL FORM. 

FOR EXAMPLE IF A MATRIX HAS 10 
COLUMNS THEN MM(1)=A(1,1) 
MM(2)=A(1,2),MM(11)=A(2» 1)» ETC. 


DIMENSION VV( 300 ) * SLL ( 300 ) » XNF( 300) » XNNUC300 ) , XNNL ( 300 ) , 
»SIGL ( 300 ) , SIGU( 300 ) , UU( 300 ) , STRNK 300 ) , STRN2( 300 ) , STMEN1 (300 
x ) , STMEN2C 300 ) , RQ( 300 ) , RMC 300) 

X, SMAX1 ( 300 ) , SMIN1 ( 300 ) , SMXMN( 2) 

x, MM(900) , RAC 32) , CMM( 32) , X2C900) , Y2C900) , BUC900 ) , BL (90 
xo ), X3( 1999 >, Y3( 1999 ),MMM( 1999 >,X4( 1999 ),Y4(1999),SD( 300) 
x,STRSl(300),STR$2(300)»DMU(300),DML(300) 

REAL KT,KT1 
READ(5,X)L 
READ(5,x)KT 
READ(5,»)AMAX 

READ(5,X)RMMIN,LEVLM,XIM,RMIN,LEVLR,XIR 
READ(5,X)EM,A,SH,SFP,BS,EFP,CS 
WRITE(6» 300)EM, SFP, BS» EFP,CS, A,SH,AMAX 
READ ( 5 , x ) ( MMM( I ) , I = 1 , L ) 

CMM(1)=RMMIN 
DO 904 L=2, LEVLM 
LL=L-1 

CMM(L)=CMM(LL)+XIM 

904 CONTINUE 
RAC 1 ) =RMIN 

DO 905 L=2> LEVLR 
LL=L-1 

905 RA(L)=RA(LL)+XIR 
R=1/SH 

IKK = 1 

DO 907 I=1,LEVLR 
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22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 


53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 


77 


DO 908 J=1 , LEVLM 
X4(IKK)=((2.XCMM(J))+RA<I))/2. 

Y4(IKK)=(2. XCMM( J ) )-X4( IKK) 

X3( IKK)=X4( IKK)XAMAX 
Y3(IKK)=Y4( IKK)XAMAX 
IKK=IKK+1 

908 CONTINUE 
907 CONTINUE 

IKK=IKK-1 
IK = 1 

DO 960 1=1, IKK 

IF(MMM(I) .EQ.O) GO TO 960 

X2(IK)=X3(I) 

Y2( IK) =Y3( I ) 

MM( IK) =MMM( I ) 

IK=IK+1 

960 CONTINUE 
IK=IK-1 
SMAX=X2( IK) 

SMIN=Y2C IK) 

DO 909 1=1, IK 

909 SLL(I)=CX2(I)-Y2(I))/2. 

NLOAD=IK 

DO 911 1=1, IK 
SMAX1(I)=X2(I) 

IF(SMAXICI) .GT.Y2CI)) GO TO 912 
SMAX1 C I )=Y2( I ) 

SMINKI )=X2(I) 

GO TO 911 

912 SMIN1(I)=Y2(I) 

911 CONTINUE 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C DETERMINING MAX, MIN LOCAL STRAIN AND STRESS 
IF(ABS(SMIN) . LT . ABS(SMAX) )THEN 
C=((SMAXxKT)xx2)/EM 
CALL NEUBER(C, R, A, V, U, EM) 

FMAX=U 

STMAX=V 

SSS=SLL(IK) 

C=C(SSSXKT)XX2)/EM 
CALL NEUBER(C,R,A,V,U,EM) 

FMIN=FMAX-(2.XU) 

STMIN=STMAX-(2.XV) 

ELSE 

C=((SMINxKT)xx2)/EM 
CALL NEUBER(C, R, A, V, U, EM) 

FMIN=U 

STMIN=V 

SSS=SLL(IK) 

C=((SSSXKT)XX2)/EM 
CALL NEUBER(C,R,A,V,U,EM) 

WRITE(6 , x)V, U 
FMAX=FMIN+(2.xU) 

STMAX=STMIN+C2.XV) 

END IF 

WRITE(6,310) FMAX, STMAX 
WRITE(6 , 320) FMIN, STMIN 
Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DETERMINING UPPER AND LOWER BOUNDS ON MEAN STRESS 
DO 1100 1=1 , NLOAD 


Table A-3 (3rd page) 


78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 


105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 
123 
’24 

*i_5 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 


SPRMl=CSMAXl(I)-SMIN>/2. 

C=((SPRM1*KT)**2)/EM 
CALL NEUBER (C,R,A,V,U,EM) 

STRN1(I)=(2.XV)+STMIN 
STRSM I ) = (2. XU)+FMIN 
SPRM2=(SMAX-SMINl(I))/2. 

C=C(SPRM2XKT)X*2)/EM 
CALL NEUBER (C,R,A,V,U,EM) 

STRN2( 1 )=STMAX-(2 . XV) 

STRS2(I)=FMAX-(2.XU) 

SL=SLL(I) 

C=CCSLXKT)XX2)/EM 

CALL NEUBER(C,R,A,V,U,EM) 

VV(I)=V 

UUCI)=U 

STMENld )=STRN1( I )-VV( I ) 

STMENZa ) =STRN2( I )+VVC I ) 

1100 CONTINUE 

DO 1310 I=l,NLOAD 
SIGUC I ) = STR51d )-UU( I ) 

SIGL ( I ) =STRS2( I )+UU( I ) 

1310 CONTINUE 

DO 1099 1=1, IK 
1099 SD(I)=2.XSLLCI) 

HRITE(6,330) 

DO 1400 I=l,NLOAD 
ER=VV( I ) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C UPPER/LOWER LIFE PREDICTION 

CALL LIFE <ER,SFP,BS,EFP,CS,EM,XNFF) 

XNFd)=XNFF 
1400 CONTINUE 

DO 1500 1=1 , NLOAD 

XNNUCI)=(XNF(I))/((1-(SIGL(I)/SFP))XXC1/BS)) 

XNNL(I)=(XNF(I))/((1-(SIGU(I)/SFP))XX(1/BS)) 

1500 CONTINUE 

DO 1600 1=1, IK 
BUd)=MMd)/XNNLd) 

BL C I ) =MM( I )/XNNU( I ) 

1600 CONTINUE 
BUU=0 . 

BLL = 0 . 

DO 1800 1=1, IK 
BUU=BU( I J+BUU 
BLL=BL(I)+BLL 
1800 CONTINUE 

BLLL=1./BLL 
BUUU=1./BUU 
DO 2001 1 = 1, IK 
DML(I)=CBL(I)/BLL)X100. 

I F C DML ( I ) . L T . 0 . 0 1 ) DML ( I ) = 0 . 

2001 CONTINUE 

DO 2200 1=1, IK 
DMU(I)=(BU(I)/BUU)X100. 

I F ( DMU C I ) . L T . 0 . 0 1 ) DMU ( I ) = 0 . 

2200 CONTINUE 

DO 1700 1=1, IK 
RQ(I)=ABS(X2( I )-Y2( I ) ) 

RMCI)=(X2(I)+Y2(I))/2. 

WRITEC6,34)RQCI),RM(I), VV(I) ,MM(I) ,SIGL(I) ,SIGU(I) , DML( I ) , DMU( I ) 
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134 1700 CONTINUE 

137 WRITE?6,350)BUUU 

138 WRITE?6,360)BLLL 

139 300 FORMAT? '1 ' ,/12X, 'ELASTIC MODUlUStKSI) 3 ' , F7 . 0//12X, 

x'FATIGUE STRENGTH COEFFICIENTtKSI) 3 • , F6 . 1//12X, 

** FATIGUE STRENGTH EXPONENT = • , F7 . A//12X, 

*' FATIGUE DUCTILITY COEFFICIENT =',F8.A//12X, 

*' FATIGUE DUCTILITY EXPONENT = ' , F7 . A//12X, 

x'CYCLIC STRENGTH COEFFICIENTtKSI) = • , F6 . 1//12X, 
x'CYCLIC STRAIN HARDENING EXPONENT =',F6.3//12X, 
x' LOAD SCALE FACTOR =',F6.3) 

140 310 FORMATt ' 1 ' ,//, 12X, 'MAX LOCAL STRESStKSI) 3 ' , F9 . 3,//, 12X, 

X'MAX LOCAL STRAIN 3 ' , F9 . 5) 

141 320 FORMAT?//, 12X, 'MIN LOCAL STRESS?KSI) 3, ,F9.3,//, 12X, 

x 'MIN LOCAL STRAIN', F9. 5) 

142 330 x ^^{JJ^2X,'N0MINAU,3X,'N0MINAL*,21X, 'LOWER', AX, -UPPER', 3X, 

X2X 3 ?0F S T0?AL' '4X, *ST*ESS',21X, ’MEAN',5X, 'MEAN', AX, 'OF TOTAL', 

x/lAX, 'RANGE', 5X, 'MEAN* , AX, 'STRAIN' , AX, 'NO OF' , 3X, 'STRESS' , 3X, 
x'STRESS',2X, 'LOWER', 5X, 'UPPER', 

X/14X, ' ?KSI ) ' , 5X, » t KSI ) ' , 2X, 'AMPLITUDE' ,2X, 'CYCLES' , 2X, ' tKSI ) ' , 
XAX, • tKSI )', 3X, • DAMAGE' , AX, 'DAMAGE* ) 

1A3 3A FORMAT ? 12X, F7 . 2, 3X, F7 . 2, 3X, F8 . 5, 2X, I5,2X,F7 . 3,2X, F7 . 3, 3X, F5 . 2, 

*5X, F5 . 2 ) 

1AA 350 FORMAT?//, 12X, 'LOWER LIFE BOUND 3 ', Ell. 3) 

1A5 360 FORMAT?//, 12X, 'UPPER LIFE BOUND 3 ', Ell .3) 

1A6 STOP 

1A7 END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXK 


1A8 

1A9 

150 

151 

152 

153 
15A 

155 

156 

157 

158 

159 

160 
161 
162 
163 
16A 

165 

166 
167 


SUBROUTINE NEUBERIC, R, A, V,U, EM) 

DIMENSION EtlOO) 

Sl = t EM/AxxR)xx( 1/? 1-R) ) 

EP1=?S1/A)XXR 

C1=SIXEP1 

IF(C.GT.Cl) GO TO 10 
Et 1 ) = ?C/EM)XX . 5 
GO TO 20 

10 E?1)=?C/A)X?R/?1+R)) 

20 DO 1 1=1,100 

P=C/EM/E?I) 

Q=tC/A/EtI))xxR 

E?I+1) 3 E?I)-?E?I)-P-Q)/?1+P/E?I)+RXQ/E?I)) 
X=E? 1+1 )/Et I ) 

IF?X. GT. 0.999. AND. X.LT. 1.001) GO TO 30 
1 CONTINUE 

30 V=Et 1+1 ) 

U=C/V 

RETURN 

END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


168 

169 

170 

171 

172 

173 
17A 


SUBROUTINE LIFE? ER, SFP, BS, EFP,CS, EM,XNFF) 
DIMENSION XNtlOO) 

XNT 3 . 5x(SFP/EM/EFP)xx? 1/tCS-BS) ) 
EPT=EFPx?2xXNT)xxCS 
IF?ER.GT.EPT)GO TO AO 
XNtl)=.5xtEMxER/SFP)xxtl/BS) 

GO TO 50 
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175 

176 

177 

178 

179 


180 

181 

182 

183 

189 

185 


90 XN<1)=.5*<ER/EFP)*«(1/CS) 

50 DO 210 K=1 , 100 

Q*SFP /EH»2» «BS 
H=EFP*2**CS 

XN ( K+l)^XN(JOr ( G*XN C K ) *«S+H*XN ( K ) JWCS-ER )✓ < OUBSKXH ( K ) *K < BS- 1 ) +H 
a*CS*XN(K)**CCS-l)) 

Z=XN(K+1)/XN(K) 

IF(Z. GT. 0.999. AND. Z.LT. 1.0010) GO TO 60 
210 CONTINUE 
60 XNFF=XN(K+1) 

RETURN 

END 

C9ENTRY 
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ELASTIC MODULUS(KSI)= 15838. 

FATIGUE STRENGTH COEFFICIENT(KSI) » 
FATIGUE STRENGTH EXPONENT = 

FATIGUE DUCTILITY COEFFICIENT 
FATIGUE DUCTILITY EXPONENT 
CYCLIC STRENGTH COEFFICIENT(KSI) * 
CYCLIC STRAIN HARDENING EXPONENT » 
LOAD SCALE FACTOR 


220.9 

-0.0763 

6.2160 

- 1.0101 

192.4 

0.076 

59.000 
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MAX LOCAL STRESS(KSI)- 124.571 
MAX LOCAL STRAIN- 0.01103 


MIN LOCAL STRESS(KSI)- -51.903 
MIN LOCAL STRAIN -0.00018 


NOMINAL 

STRESS 

RANGE 

NOMINAL 

STRESS 

MEAN 

STRAIN 

NO OF 

LONER 

MEAN 

STRESS 

UPPER 

MEAN 

STRESS 

PERCENT 
OF TOTAL 
LONER 

PERCENT 
OF TOTAL 
UPPER 

(KSI) 

(KSI) 

AMPLITUDE 

CYCLES 

(KSI) 

(KSI) 

DAMAGE 

DAMAGE 

2.34 

30.48 

0.00019 

1 

53.778 

54.296 

0.00 

0.00 

2.34 

35.40 

0.00019 

2 

65.578 

66.093 

0.00 

0.00 

2.34 

37.74 

0.00019 

16 

71.478 

71.990 

0.00 

0.00 

2.34 

40.12 

0.00019 

2 

77.378 

77.885 

0.00 

0.00 

4.72 

33.04 

0.00037 

1 

59.678 

60.194 

0.00 

0.00 

4.72 

37.74 

0.00037 

4 

71.478 

71.988 

0.00 

0.00 

4.72 

40.12 

0.00037 

2 

77.378 

77.881 

0.00 

0.00 

7.08 

37.74 

0.00056 

2 

71.478 

71.985 

0.00 

0.00 

11.80 

35.40 

0.00093 

1 

65.578 

66.085 

0.00 

0.00 

14.52 

37.74 

0.00130 

2 

71.478 

71.960 

0.00 

0.00 

21.24 

30.48 

0.00168 

1 

53.778 

54.285 

0.00 

0.00 

21.24 

40.12 

0.00168 

1 

77.378 

77.793 

0.00 

0.00 

23.40 

25.94 

0.00186 

12 

41.979 

42.492 

0.00 

0.00 

23.40 

30.68 

0.00186 

43 

53.778 

54.281 

0.00 

0.00 

23.40 

33.04 

0.00186 

30 

59.678 

60.169 

0.00 

0.00 

23.40 

35.40 

0.00186 

1742 

65.578 

66.049 

0.11 

0.11 

23.40 

37.74 

0.00186 

1348 

71.478 

71.916 

0.14 

0.14 

23.40 

40.12 

0.00186 

27 

77 .378 

77.765 

0.00 

0.00 

23.40 

42.48 

0.00186 

223 

83.278 

83.584 

0.07 

0.07 

25.94 

30.68 

0.00205 

1 

53.778 

54.276 

0.00 

0.00 

25.94 

35.40 

0.00205 

1 

65.578 

66.034 

0.00 

0.00 

25.94 

40.12 

0.00205 

1 

77.378 

77.729 

0.00 

0.00 

28.32 

30.68 

0.00224 

14 

53.779 

54.269 

0.00 

0.00 

28.32 

33.04 

0.00224 

5 

59.678 

60.149 

0.00 

0.00 

28.32 

35.40 

0.00224 

453 

65.578 

66.016 

0.30 

0.31 

28.32 

37.76 

0.00224 

2613 

71.478 

71.865 

2.89 

2.92 

28.32 

40.12 

0.00224 

155 

77.378 

77.684 

0.29 

0.29 

28.32 

42.48 

0.00224 

12 

83.278 

83.461 

0.04 

0.04 

30.48 

35.40 

0.00242 

5 

65.578 

65.993 

0.00 

0.00 

30.48 

37.76 

0.00242 

20 

71.478 

71.829 

0.06 

0.06 

30.48 

40.12 

0.00242 

1 

77.378 

77.629 

0.00 

0.00 

33.04 

30.68 

0.00261 

24 

53.779 

54.249 

0.05 

0.05 

33.04 

33.04 

0.00261 

6 

59.679 

60.116 

0.02 

0.02 

33.04 

35.40 

0.00261 

65 

65.578 

65.965 

0.33 

0.33 

33.04 

37.76 

0.00261 

7785 

71.478 

71.784 

64.88 

65.23 

33.04 

40.12 

0.00261 

460 

77.378 

77.561 

6.50 

6.47 

33.04 

42.48 

0.00261 

4 

83.278 

83.278 

0.10 

0.10 

35.40 

33.04 

0.00279 

1 

59.679 

60.093 

0.00 

0.00 

35.40 

37.76 

0.00279 

15 

71.478 

71.729 

0.31 

0.31 

37.74 

30.68 

0.00298 

6 

53.780 

54.216 

0.07 

0.07 

37.74 

33.04 

0.00298 

8 

59.679 

60.065 

0.14 

0.14 

37.74 

35.40 

0.00298 

30 

65.578 

65.884 

0.87 

0.87 

37.74 

37.76 

0.00298 

6 

71.478 

71.661 

0.29 

0.29 

37.74 

40.12 

0.00298 

24 

77.378 

77.378 

1.95 

1.91 

40.12 

30.68 

0.00317 

1 

53.781 

54.194 

0.02 

0.02 
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40.12 

33.04 

42.48 

30.68 

42.48 

33.04 

42.48 

35.40 

44.84 

30.68 

44.84 

33.04 

44.84 

35.40 

47.20 

30.68 

47.20 

33.04 

47.20 

35.40 

51.92 

33.04 

70.80 

23.60 

LOHER LIFE 

BOUND* 

UPPER LIFE 

BOUND* 


0.00317 

1 

0.00335 

10 

0.00335 

4 

0.00335 

28 

0.00354 

1 

0.00354 

1 

0.00354 

1 

0.00373 

1 

0.00373 

1 

0.00373 

5 

0.00410 

1 

0.00560 

1 

0.851E 03 


0.869E 03 



59.690 60.029 

53.782 56.165 

59.680 59.984 

65.579 65.762 

53.784 54.129 

59.681 59.930 

65.579 65.679 

53.785 54.085 

59.681 59.863 

65.579 65.579 

59.682 59.682 

36.337 36.337 


0.04 0.04 

0.52 0.52 

0.33 0.33 

3.78 3.76 

0.11 0.11 

0.17 0.17 

0.27 0.27 

0.21 0.21 

0.33 0.33 

2.69 2.63 

1.15 1.13 

10.96 10.73 
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COMPUTER PROGRAM FOR RAIN-FLOW CYCLE COUNTING ANALYSIS 
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( WP : comp uter- 2) 



ABSTRACT 


A computer program (RAINF) for rain-flow cycle counting analysis is 
provided. The program can take a lengthy load history and reduce it to 
the compact form of a matrix giving combinations of range and mean or 
peak and valley values. This information can be used for fatigue analy- 
sis. The input values are defined and five examples using different 
options of the program are provided. 
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INTRODUCTION 


Rain-flow cycle counting is a method that exists for reducing 
service load history records to a compact description so that the infor- 
mation can be used in analysis for fatigue. The compact description is 
in the form of a matrix giving combinations of range and mean, or peak 
and valley values. This method is widely accepted as the most accurate 
cycle counting method for predicting fatigue life based on a cumulative 
damage type of approach. The rain-flow method cannot be misled by any 
synthetic load sequences and will always count the cycles correctly, 
based on the fact that closed hysteresis loops are most representative 
of a fatigue damage event. 

Figure 8-1 illustrates rain-flow cycle counting. As shown in this 
example, the rain-flow method has the important characteristic that it 
counts the major load excursions as cycles, while also counting the 
minor events. This feature allows it to realistically handle real 
service loading where there are low level vibratory loads, etc., 
superimposed on major cycles associated with the usage of the machine, 
vehicle or structure, such as ground-air-ground cycles in aircraft. 

The computer program explained in this manual is the same program 
that was previously used in Reference B-l. 
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Figure B-l. Illustration of rain-flow cycle counting from the ASTM 
Standards [B-21. 










PROGRAM LOGIC 


The following logic is used consistent with the ASTM standard 
[B-2]: Let x denote the absolute value of the range under consideration, 
and y previous absolute range adjacent to x. 

Step 1: Determine the maximum absolute value in the history. (Note 

that this value can be either a peak or a valley.) 

Step 2: Arrange the history to start with the maximum absolute value. 

Move all peaks and valleys which occur prior to the maximum to the end 
as Illustrated in Figure B-la. 

Step 3: Read the next value. If out of data, go to step 9. 

Step 4: Three points are needed to define x and y. If there are less 

than three points, go back to step 3. Define x and y using the three 
most recent peaks and valleys that have not been discarded. 

Step 5: Compare the two ranges, namely x and y. 

a) If x < y, go to step 3. b) If x > y, go to step 6. 

Step 6: If a rain-flow filtered history is not desired, go to step 8. 

Step 7: If y < filter level specified in the program Input, then 
discard the peak and valley of the range y in the array in memory which 
is the original history of step 2. 

Step 8: Count range y as one cycle, determine the mean value of the 

peak and valley of y, discard the peak and valley of y in the array set 

up in step 3, and go to step 4. 

Step 9: Stop. 
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DEFINITION OF INPUT DATA 


Line No. of 

READ Statement Variable Name Explanation 


3 


OPTION 


= 1 List filtered history 
as peak/valley sequence; 
also print range/mean 
matrix of rain-flow 
cycles for original 
history. 

= 2 List range, mean, min, 
and a max of rain-flow 
cycles not in matrix form. 

= 3 Print range/mean matrix of 
rain-flow cycles. 

= 4 Print max/min 32x32 matrix of 
rain-flow cycles. 


4 

FL 

Filter level as a range 
value 

5 

NN 

Number of peak/valley points in 
history 

8 

XIM 

Constant increment between mean 
values in the range/mean matrix 


XIR 

Constant increment between range 
values in the range/mean matrix 

9 

P( ) 

Input load history as peaks and 
valleys in sequence 


Comment 


Integer 


The history is con- 
verted to a min 
value of 1 and max 
value of 32. 

Required only for 
OPTION = 1. 


Required only for 
OPTION = 1 or 3. 

Required only for 
OPTION = 1 or 3. 

History must start 
and end with the 
same value. 
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EXAMPLE 1 


The history of Figure B-l is used for this rain-flow cycle counting 
example. Option 2 of the program is used; therefore, the result is 
shown as a list of range, mean, minimum, and maximum values. Table B-l 
shows the input for this example. The entire program listing and 
program output for this example are attached as Table B-2. 
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Table B-l. Input for Example 1 



Table B-2. Program Listing and Output for Example 1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


C$ JOB 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


40 


NPV1336,P=100 

RAIN-FLOW COUNTING PROGRAM (RAINF) 

NOTE THAT THE HISTORY MUST START AND END WITH THE SAME VALUE. 


INPUT 

DATA LINE 1 


DATA LINE 2 
DATA LINE 3 
DATA LINE 4 


DATA LINE 5 


OPTIONS LIST FILTER HISTORY AS PEAK/VALLEY 
SEQUENCE. ALSO PRINT RANGE/MEAN 
MATRIX OF RAINFLOW CYCLES FOR 
ORIGINAL HISTORY. 

=2 LIST RANGE, MEAN, MIN, AND MAX OF RAIN- 
FLOW CYCLES NOT IN MATRIX FORM. 

=3 PRINT RANGE/MEAN MATRIX OF RAINFLOW 
CYCLES 

=4 PRINT MAX/MIN MATRIX OF RAINFLOH 
CYCLES . 

FL=FILTER VALUE AS A RANGE 
NN=NUMBER OF PEAK/VALLEY POINTS IN HISTORY 
XIM=CONSTANT INCREMENT BETWEEN MEAN VALUES 
IN THE RANGE/MEAN MATRIX. _ 

XIR=CONSTANT INCREMENT BETWEEN RANGE VALUES 
IN THE RANGE/MEAN MATRIX. 

NOTE THAT XIM,AND XIR REQUIRED FOR OPTIONS OR 3 

P( )=INPUT LOAD HISTORY AS PEAKS AND VALLEYS 
IN SEQUENCE. 


REAL P< 10000) ,PE( 10000) ,PP( 10000), PC(10000>,PCC(10000),MM«4>, 
*RA< 64), PI ( 10000 ),MEAN( 5005), R( 5005) 

INTEGER M(64,64),SUM(64),SUMM,MI<32, 32), OPTION 
READ(5,X) OPTION 
I F( OPTION . EQ . 1 )READ(5, *)FL 
READ(5,*)NN 

I F( OPTION. EQ.4)G0 TO 40 
IF(0PTI0N.EQ.2)G0 TO 40 
READ(5,X)XIR,XIM 
READ(5,*HP(I),I = i,NN) 


N=NN 


C 

C 

C DETERMINATION OF LARGEST PEAK OR VALLEY 
LCOUNT=l 
DO 100 1=1, N 
PECI )=P(I) 

100 CONTINUE 

PMAX=ABS(P(1)) 

DO 200 1=2, N 

IF(PMAX.GE.PEd) ) GO TO 200 
PMAX=ABS(PE( I ) ) 

LCOUNT=I 
200 CONTINUE 

I F( OPT I ON . EQ . 4)THEN 
SMAX=P( 1 ) 

SMIN=P( 1 ) 

DO 301 1=2, NN 
IF(P(I).GT. SMAX)SMAX=P( I ) 

IFCPCI) . LT . SMIN)SMIN=P( I ) 

301 CONTINUE 
CF1=SMIN 
CF2=SMAX 
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Table B-2 (2nd page) 


30 


CF3=SMAX-SMIN 

31 


END IF 


C 



C 



C 

ARRANGE THE PEAK OR VALLEY 

32 


JK=LC0UNT+1 

33 


J=N-JK+1 

34 


KKK=LCOUNT 

35 


DO 300 1=1, J 

36 


PP( I )=P(KKK) 

37 


KKK=KKK+1 

38 

300 

CONTINUE 


C 

JJJ=N-LC0UNT-1 

39 


J=J+1 


C 

DO 130 1=1, JJJ 

40 


DO 350 1=1 , LCOUNT 

41 


PP(J)=P(I) 

42 


J=J+1 

43 

350 

CONTINUE 

44 


DO 500 1=1, NN 

45 


PC(I)=PP<I) 

46 

500 

CONTINUE 

4 T 


NNN=N+1 

48 


I F( OPT ION . EQ .2)WRITE(6 , 210) 


C 



C 



C 

FINDING THE CYCLE 

49 


AA=3 . 1422 

50 


DO 194 1=1,32 

51 


DO 195 J=1 , 32 

52 


MICI, J)=0 

53 

195 

CONTINUE 

54 

194 

CONTINUE 

55 


1 = 0 

56 


K = 1 

57 


IF(0PTI0N.EQ.2)WRITE(6,107) 

58 


J=1 

59 

2 

1 = 1+1 

60 


IFd . LT.3) GO TO 2 

61 


J = J+1 

62 


IF(I.E9.NNN) GO TO 400 

63 

50 

IFCPP(J).EQ.AA) THEN 

64 


J = J-1 

65 


GO TO 50 

66 


END IF 

67 


JM1=J-1 

68 

60 

IF(PPCJMl) . EQ.AA) THEN 

69 


JM1=JM1-1 

70 


GO TO 60 

71 


END IF 

72 

70 

IF(I.GT.NNN) GO TO 400 

73 


X=ABS(PP(I )-PP( J) ) 

74 


Y=ABS(PP(J)-PP(JM1)) 

75 


XX=(PP(J)+PP(JMl))/2. 

76 

5 

IF(X.GE.Y)THEN 

77 


I F< OPTION . NE . 1 )G0 T01600 

78 


IF(Y.LE.FL) THEN 

79 


PCCJ)=AA 

80 


PC< JM1 )=AA 

81 


END IF 
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Table B-Z (3rd page) 


82 

1600 

I F< OPTION . EQ .2)60 TO 41 

83 


I F( OPT I ON . EQ . 4)THEN 

8$ 


PI(J) = U32.*(-CF1+PP(J)))+CF2-PP(J))/CF3 

85 


PI(JM1)=((32. X(-CF1+PP( JM1) ))+CF2-PP( JM1))/CF3 

86 


PMAX=PI(J) 

87 


PMIN=PI( JM1) 

88 


IF(PMAX.LT.PMIN)THEN 

89 


PMAX=PI( JM1) 

90 


PMIN=PI(J) 

91 


END IF 

92 


END IF 

93 


I F( OPT ION . EQ . 4JTHEN 

94 


El J=PMAX+ . 5 

95 


EJI=PMIN+ . 5 

96 


IJ=INT(EIJ) 

97 


JI=INT(EJI) 

98 


MI(IJ,JI)=MI(IJ,JI)+1 

99 


END IF 

100 


I F( OPT I ON . EQ . 4)60 TO 42 

101 


R( K) =Y 

102 


MEAN(K)=XX 

103 


K=K+1 

104 


GO TO 42 

105 

41 

PMAX=PP( J ) 

106 


PMIN=PP( JM1 ) 

107 


IF(PMAX. IT . PMIN)THEN 

108 


PMAX=PP( JM1 ) 

109 


PMIN=PP(J) 

110 


END IF 

111 


WRITE(6 , 108 >Y, XX, PMAX, PMIN 

112 

42 

PP( J)=AA 

113 


PP( JM1 ) =AA 

114 


J=J-1 

115 

3 

J=J-1 

116 


IF(J.LT.l) GO TO 11 

117 


IF(PP( J ) . EQ . AA) GO TO 3 

118 


IF( I . EQ . 4) THEN 

119 


1=5 

120 


J = I-1 

121 


END IF 

122 


JM1=J-1 

123 


IF(JM1.LT.1)THEN 

124 


JM1=J 

125 


J = I 

126 


1 = 1 + 1 

127 


GO TO 70 

128 


END IF 

129 

6 

IF(PP( JM1) .EQ.AA) GO TO 4 

130 


X=ABS(PP(I)-PP(J)) 

131 


Y=ABS(PP(J)-PP(JM1)) 

132 


XX=( PP( J )+PP( JM1 ) )/2 . 

133 


GO TO 5 

134 

4 

JM1=JM1-1 

135 


IF( JM1 . LT . 1 )THEN 

136 


JM1 = J 

137 


J = I 

138 


1 = 1 + 1 

139 


END IF 

140 


GO TO 6 

141 

11 

1 = 1+2 
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Table B-2 (4th page) 


142 

143 

144 

145 

146 

147 

148 

149 


150 

151 

152 

153 

154 

155 


158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 
189 


190 

191 

192 

193 

194 

195 


J = I-1 
JM1 = J-l 
GO TO 70 
ELSE 

XX=(PPCJ)+PP(JMl))/2. 

J=I-1 
END IF 
GO TO 2 
C 
C 

c 

400 IFC0PTI0N.EQ.2) GO TO 999 
IF( OPTION. EQ .4)00 TO 998 
K=K~l 
RMAX=R(1) 

RMIN=R(1) 

RMMAX=MEAN(1) 

RMMIN=MEAN( 1 ) 

DO 1800 1=2, K 
IF(R<I).GT.RMAX)RMAX=R(I) 
IF(R(I).LT. RMIN)RMIN=R( I) 

I F( MEAN ( I ) . GT . RMMAX ) RMMAX=MEAN( I ) 
IFIMEANTI} . LT. RMMIN)RMMIN=MEAN(I) 
1800 CONTINUE 

DIFR=RMAX-RMIN 

DIFM=RMMAX-RMMIN 

ER=DIFR/XIR 

EM=DIFM/XIM 

ER=ER+2 

EM=EM+2 

LEVLM=INT(EM) 

LEVLR=INT(ER) 

DO 192 L=l, LEVLR 
DO 193 LL=1, LEVLM 
M(L , LL )=0 
193 CONTINUE 
192 CONTINUE 

YA=RMIN-XIR 

XA=RMMIN-XIM 

WRITE(6, 111 )RMIN,RMAX,RMMIN, RMMAX 
XB= . 50 
YB= . 50 

DO 1900 1=1, K 
EI=((R(I)-YA)/XIR)+YB 
EJ=((MEAN(I)-XA)/XIM)+XB 
II=INT(EI) 

JJ=INT(EJ) 

IFCII.EQ. 0)11=1 
IF(JJ.EQ.0)JJ=1 
M( II , J J ) =M( II , J J )+l 
1900 CONTINUE 
C 
C 

C FILTERING PROCESS 

I F( OPT I ON . NE. 1 )G0 TO 1102 
KN = 1 

DO 1000 11=1, NN 

IF(PCCII) .EO.AA) GO TO 1000 

PCC<KN)=PCCII) 

KN=KN+1 
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196 

1000 

197 

198 

199 

112 

200 

201 

113 

202 

203 

1103 

204 

205 

1001 

206 

207 

C 

C 

C 

998 

208 

209 

210 
211 
212 

213 

214 

215 

216 

202 

217 

201 

218 

1102 

219 

220 
221 
222 
223 

900 

224 

225 

226 
227 

1100 

228 

229 

99 

230 

• 

231 

232 

233 

234 

98 

235 

236 

1153 

237 

999 

238 

1151 

239 

240 

1152 

241 

242 

243 

996 

244 

245 

1154 

246 

247 

604 

248 

605 

249 

250 

1150 

251 

600 


CONTINUE 
KN=KN-1 
MRITEC6 112) 

FORMAT< '1 ' ,//15X, ’FILTER HISTIRY-PEAK/VALLEY SEQUENCE') 
WRITEC6 , 113)FL 

FORMATL//15X, 'FILTER LEVEL*' , F7 , 3) 

F0RMAT(//15x! 'NUMBER OF POINTS IN FILTER HISTORY* M5,//) 
WRITE!6,1001)!PCC!I),I=1,KN) 

FORMAT! 1X,8! 2X> F6 . 1 ) ) 


MATRIX PREPRATION 

GO TO 1102 

LEVLM=32 

LEVLR*32 

XIR=1 

XIM=1 

RMMIN=1 

RMIN=1 

DO 201 1*1,32 
DO 202 J=1 , 32 
MCI, J)*MI(I, J) 

CONTINUE 

CONTINUE 

I F( OPT I ON . EQ . 2)G0 TO 999 
MM( 1 ) =RMMIN 
DO 900 L=2, LEVLM 
L L = L - 1 

MMC L ) =MM! LL )+XIM 

CONTINUE 

RA( 1 )=RMIN 

DO 1100 L=2, LEVLR 

LL=L-1 

RA(L)=RA(LL)+XIR 
1 = 0 
1 = 1+1 

IF(I .GT. LEVLR) GO TO 1153 
SUM(I)=0. 

DO 98 J=l, LEVLM 

SUM! I ) =SUMC I )+M( I » J ) 

CONTINUE 

GO TO 99 

CONTINUE 

I F( OPT I ON . EQ . 2)G0 TO 997 


L = 1 


LB = 8 


I F( OPTION . EQ . 1 )G0 TO 996 
HRITE!6 , 116 ) 

GO TO 1154 
HRITE!6 , 114) 

WRITE! 6 , 1 15) FL 
I F( OPT I ON . EQ . 4)G0 TO 604 
GO TO 1150 
WRITE! 6 , 605) 

F0RMATC77X, 'TOTAL' , /IX, 'PEAK 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


/XXXXXXXXXXXXXXXXXXXXXXXXXXXXXVAL L EY 
' , 4X, 'CYCLES' ,/) 


GO TO 2100 

FORMAT! 77X? 'TOTAL ' ,/lX, 'RANGE /xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxMEA 
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252 

2100 

253 


254 


255 


256 

1300 

257 


258 


259 


260 


261 


262 

1400 

263 


264 

1500 

265 


266 

101 

267 

102 

268 

103 

269 

104 

270 

997 

271 

107 

272 

108 

273 

111 

274 

116 

275 

114 

276 

115 

277 

210 

278 


279 



xnxxxxxxkxxxxxxxxxxxxxxxxxkxxxxx • , 3X. 'CYCLES' ./) 
WRITEC6 , 101 ) (MMCLL ) , LL=L , LB) 

DO 1300 1=1 , LEVLR 
WRITE(6,102)RA(I),SUM<I) 

WRITE(6 , 103HM(I,J),J=L,LB) 

CONTINUE 

IF( LB . EQ. LEVLN)GO TO 1400 

L =L+8 

LB=LB+8 


IFCLB.GT. LEVLM)LB=LEVLM 

GO TO 1152 

SUMM=0 

DO 1500 1=1, LEVLR 
SUMM=SUMM+SUM(I) 

WRITEC6 , 104)SUMM 
FORMAT< 12X,8(F6 . 1 , 2X) ) 

FORMAT (2X, F6.1,69X,I4) 

FORMAT ( '+',11X,8(I6,2X)) 

FORMAT(//,5X, 'TOTAL NO OF CYCLES* ', 3X, 15) 

CONTINUE 

F0RMAT(15X, 'RANGE' , 15X, 'MEAN* , 15X, ' MAX',15X,» MIN') 

FORMAT ( 14X, F7 . 3, 12X, F7 . 3, 13X, F7 . 3, 13X, F7* 3) 

FORMAT C ' 1 ' ,//15X, 'MIN RANGE* ' , F8 . 3,//15X, 'MAX RANGE=',F8.3, 

X//15X, 'MIN MEAN=',F8.3,//15X,'MAX MEAN*' ,F8.3) 

FORMATL ' 1 * ,//35X, 'RA INFLOW CYCLES ') 

FORMAT('1',//20X, 'RAINFLOW CYCLES FOR ORIGINAL HISTORY') 
F0RMATC/5X, *NO RANGE LESS THAN OR EQUAL TO FILTER LEVEL*', F7. 3 
x,3X, 'OCCUR IN FILTER HISTORY*) 

xing M method'^) 10X ' ,RANGES C0UNTEDL AS cycles ® y RAINFLOH CYCLE COUNT 

STOP 

END 


C$ ENTRY 
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RANGES COUNTED AS CYCLES BY RAINFLOH CYCLE COUNTING METHOD. 


RANGE 

4.000 

3.000 

7.000 

9.000 


MEAN 

1.000 

- 0.500 

0.500 

0.500 


MAX 

3.000 

1.000 

4.000 

5.000 


MIN 

- 1.000 

- 2.000 

- 3.000 

- 4.000 
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EXAMPLE 2 


A history containing 41 peak/valley points is used with Option 2. 
Tables B-3 and B-4 show the input and output of this program, 
respectively. 
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Table B-3. Input for Example 2 


2 

41 

-8. 6, 1.11, -6 .2,-1 .71,-2.95,2. 10, -6. 6, 5. 43, -2. 58, -.4, -2. 9, 6'. 51,-1. 12, 
-1.65,1.88,-4.47,2.01,-2.5,7.3,-1 .08, 5. 36, -5. 87, -.18, -3. 5, .85,-1 89, 

5.3, .1,4. 45, -3. 82, 1.89, -2. 04, 3. 9, -4. 97, 7. 78, -3. 99, 4. 61, -5! 52, 3.1; -l! 


.97, 

6 
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Table B-4. Output for Example 2 


RANGES COUNTED AS 
RANGE 
1,240 
7,310 
8.700 
2.180 
8.330 
2.090 
3.530 
4.510 
10.980 
6.440 
3.320 
2.740 
4.350 
3.930 
7.720 
10.270 
13,170 
8.600 
8.620 
16.380 


CYCLES BY RAINFLOW 
MEAN 

- 2.330 
- 2.545 
- 2.250 
- 1 . 490 
1.265 
- 0.075 
0.115 
- 0.245 
1.020 
2.140 
- 1.840 
- 0.520 
2.275 
- 0.075 
0 . 040 
0.165 
0.715 
0.310 
- 1.210 
- 0.410 


CYCLE COUNTING METHOD. 
MAX 
- 1.710 
1.110 
2.100 
- 0.400 
5.430 
0.970 
1.880 
2.010 
6.510 
5.360 
- 0.180 
0.850 
4.450 
1.890 
3.900 

5.300 

7.300 
4.610 
3.100 
7.780 


MIN 

- 2.950 

- 6.200 

- 6.600 

- 2.580 

- 2.900 

- 1.120 

- 1.650 

- 2.500 

- 4.470 

- 1.080 

- 3.500 

- 1.890 

0.100 

- 2.040 

- 3.820 

- 4.970 

- 5.870 

- 3.990 

- 5.520 

- 8.600 
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EXAMPLE 3 


A history containing 1709 peak/valley points, called the SAE Trans- 
mission History, is used with Option 3. Therefore, the results are 
given in the form of a compact matrix containing range and mean values. 
The results are in agreement with the published values from Reference B- 
3. Tables B-5 and B-6 show the input and output for this example, 
respectively. 

Note that the size of the matrix is determined by the XIM and XIR 
increment values, for mean and range, respectively, which are given to 
the program as input data. The program determines the smallest range 
and the lowest mean value. The increments are then added by the program 
a sufficient number of times to include the largest range and highest 
mean. Hence, the program automatically makes the matrix large enough to 
include the largest range and highest mean. In this particular case, 
the matrix size resulting from XIM = XIR = 50 is 27x20. Since the 
smallest range is 200 

200 + (27 - 1)50 = 1500 

which is sufficiently large to include the largest range of 1494. 
Similarly, since the lowest mean value is -160, 

-160 + (20 - 1)50 = 790 

which is sufficiently large to include the highest mean of 753. 

Note that the number of cycles counted also appears in the 
output. This is always related to the number of peak/valley points in 
the history, NN, as follows: 


No. of cycles = 


NN - 1 


2 
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NN is always odd, since the last peak or valley, P( ), is a repetition 
of the first one, in this particular case 

No. of cycles = — = 854 
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Table B-5. Input for Example 3 


3 

1709 
50 . 50 . 


0 

513 

292 

562 

267 

585 

314 

524 

299 

513 

299 

569 

272 

535 

298 

523 

291 

680 

444 

683 

396 

678 

374 

701 

276 

607 

331 

542 

339 

572 

274 

500 

285 

565 

274 

491 

237 

496 

-76 

142 

-276 

71 

-309 

19 

-249 

137 

-93 

658 

275 

534 

240 

612 

360 

690 

315 

689 

306 

698 

417 

666 

392 

699 

374 

701 

368 

659 

345 

650 

403 

610 

393 

637 

400 

618 

285 

630 

385 

616 

374 

703 

-174 

136 

-331 

334 

98 

510 

151 

542 

84 

642 

398 

651 

-197 

44 

-197 

125 

-121 

228 

-17 

631 

387 

657 

408 

623 

218 

553 

291 

640 

408 

693 

453 

657 

356 

642 

400 

698 

330 

571 

337 

689 

395 

647 

377 

594 

227 

631 

299 

682 

411 

663 

404 

629 

306 

631 

306 

614 

409 

648 

414 

683 

374 

575 

324 

691 

379 

585 

368 

621 

371 

610 

373 

582 

-8 

240 

-227 

248 

-205 

285 

-89 

470 

248 

460 

136 

519 

317 

517 

278 

503 

179 

615 

324 

559 

219 

679 

421 

683 

-254 

84 

-222 

243 

-19 

495 

27 0 

629 

291 

56 0 

351 

657 

243 

608 

373 

663 

269 

653 

156 

663 

443 

690 

249 

612 

229 

602 

339 

738 

231 

498 

237 

614 

-167 

39 

-261 

437 

192 

503 

103 

682 

201 

589 

226 

480 

135 

561 

347 

548 

149 

494 

260 

661 

158 

467 

-184 

161 

-105 

555 

248 

449 

132 

712 

276 

529 

301 

524 

309 

564 

195 

467 

35 

378 

-294 

188 

-308 

245 

-211 

607 

41 

526 

205 

506 

218 

513 

249 

603 

261 

712 

506 

727 

194 

400 

-244 

173 

-276 

610 

405 

615 

360 

642 

403 

641 

353 

556 

340 

655 

433 

669 

422 

662 

457 

796 

403 

671 

78 

549 

285 

616 

376 

690 

18 

731 

108 

578 

282 

629 

170 

748 

181 

578 

342 

604 

261 

585 

330 

694 

349 

614 

313 

521 

290 

533 

266 

565 

291 

645 

-228 

1 

-286 

157 

-87 

725 

183 

714 

393 

621 

194 

426 

152 

561 

235 

551 

213 

691 

136 

487 

164 

549 

244 

598 

162 

507 

249 

716 

339 

733 

194 

566 

-191 

120 

-179 

108 

-119 

181 

-74 

722 

422 

639 

303 

671 

369 

715 

350 

605 

153 

765 

409 

762 

355 

712 

398 

718 

362 

679 

448 

718 

446 

698 

342 

605 

188 

498 

276 

537 

283 

522 

222 

546 

345 

570 

191 

500 

280 

658 

269 

685 

341 

551 

344 

618 

303 

612 

285 

666 

228 

534 

303 

546 

-110 

110 

-362 

88 

-199 

173 

-130 

612 

333 

680 

341 

603 

283 

619 

30 

406 

206 

470 

260 

717 

218 

487 

190 

573 

270 

476 

265 

607 

310 

784 

573 

822 

496 

825 

508 

887 

491 

769 

486 

716 

342 

594 

365 

578 

254 

593 

355 

561 

324 

573 

153 

696 

443 

651 

372 

653 

-256 

39 

-173 

620 

213 

560 

217 

559 

274 

491 

211 

571 

266 

640 

169 

572 

223 

496 

221 

573 

264 

7 58 

358 

743 

501 

777 

-161 

57 

-146 

210 

-13 

216 

-15 

373 

133 

371 

151 

624 

313 

628 

405 

698 

-202 

146 

-92 

483 

141 

551 

229 

551 

207 

533 

271 

554 

218 

532 

210 

630 

179 

489 

266 

618 

408 

650 

328 

599 

-211 

158 

-204 

68 

-179 

256 

-108 

114 

-114 

476 

228 

494 

178 

658 

218 

612 

245 

652 

416 

658 

355 

678 

428 

722 

437 

747 

542 

776 

-201 

405 

188 

707 

408 

640 

422 

648 

415 

727 

393 

610 

387 

722 

417 

720 

476 

688 

484 

691 

360 

570 

-350 

528 

216 

679 

47 5 

689 

274 

562 

286 

653 

233 

453 

253 

573 

282 

582 

356 

603 

234 

871 

494 

771 

511 

773 

481 

727 

491 

836 

566 

782 

398 

659 

329 

685 

387 

629 

218 

604 

329 

628 

-248 

274 

-206 

710 

463 

663 

365 

695 

406 

803 

437 

739 

449 

679 

452 

695 

454 

703 

26 0 

602 

-330 

34 

-350 

188 

-180 

309 

-7 

648 

417 

647 

439 

683 

438 

699 

451 

712 

382 

698 

296 

625 

385 

680 

405 

718 

457 

673 

372 

744 

414 

618 

82 

597 

191 

570 

211 

522 

212 

517 

211 

495 

232 

439 

235 

616 

256 

512 

27 0 

575 

314 

791 

369 

787 

502 

769 

468 

785 

444 

798 

501 

732 

522 

723 

462 

680 

-264 

68 

-154 

188 

-14 

242 

13 

582 

291 

618 

227 

537 

285 

555 

224 

545 

329 

588 

374 

647 

242 

561 

353 

694 

433 

753 

-237 

313 

-147 

226 

-163 

705 

329 

636 

376 

680 

383 

641 

404 

712 

367 

732 

290 

610 

318 

629 

315 

588 

353 

57 7 

265 

612 

368 

723 

-183 

267 

-208 

111 

-153 

473 

271 

576 

365 

618 

395 

605 

-283 

235 

-172 

687 

242 

521 

299 

532 
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Table B-5 (2nd page) 


283 

596 

306 

629 

392 

639 

326 

409 

640 

272 

496 

175 

470 

251 

322 

529 

265 

586 

353 

555 

-215 

164 

517 

126 

560 

165 

726 

408 

378 

744 

309 

635 

417 

677 

373 

197 

548 

181 

523 

296 

582 

113 

341 

607 

317 

604 

351 

567 

141 

-183 

199 

-28 

425 

167 

620 

400 

365 

812 

395 

695 

416 

758 

469 

27 

637 

212 

566 

275 

526 

325 

281 

533 

-204 

199 

-185 

136 

-86 

253 

506 

275 

511 

108 

516 

169 

-259 

-22 

-242 

244 

1 

222 

-17 

358 

567 

356 

599 

361 

607 

394 

299 

583 

190 

506 

235 

570 

299 

270 

619 

228 

583 

371 

593 

373 

293 

717 

417 

701 

320 

709 

385 

285 

533 

310 

535 

319 

555 

326 

102 

490 

263 

534 

45 

693 

398 


189 

489 

-367 

207 

-276 

-34 

-263 

-240 

-18 

-219 

181 

-470 

-19 

-293 

1 

566 

117 

516 

148 

415 

-406 

178 

582 

331 

835 

490 

748 

355 

419 

722 

416 

769 

442 

775 

501 

415 

683 

341 

769 

485 

738 

421 

356 

675 

426 

705 

449 

688 

479 

-320 

438 

222 

639 

251 

647 

431 

342 

705 

374 

684 

427 

698 

368 

322 

624 

368 

659 

392 

748 

355 

315 

693 

330 

607 

226 

744 

443 

258 

588 

141 

457 

221 

428 

-164 

212 

733 

318 

629 

357 

752 

458 

277 

636 

358 

562 

229 

770 

535 

-233 

288 

-100 

265 

-120 

703 

14 

164 

593 

319 

632 

398 

623 

250 

-141 

140 

-129 

544 

303 

644 

358 

447 

813 

506 

817 

537 

776 

524 

471 

789 

487 

706 

452 

669 

192 

89 

576 

362 

789 

191 

449 

151 

610 

829 

530 

734 

-318 

728 

83 

609 

876 

630 

888 

527 

797 

-239 

271 

578 

324 

675 

221 

632 

362 

-140 

510 

251 

661 

366 

695 

161 

-247 

217 

-301 

149 

-207 

157 

-74 

153 

492 

270 

607 

240 

645 

202 

-138 

508 

304 

603 

67 

626 

342 

304 

669 

337 

605 

367 

567 

280 

146 

506 

294 

577 

-11 

222 

-341 

165 

398 

176 

822 

561 

801 

553 

• 52 

834 

352 

642 

285 

587 

-222 

234 

626 

226 

480 

35 

594 

189 

-229 

336 

87 

798 

529 

753 

51 

394 

830 

464 

732 

330 

800 

452 

462 

782 

116 

653 

-223 

54 

-223 

251 

658 

192 

640 

323 

578 

281 

296 

645 

353 

664 

318 

575 

178 

412 

615 

344 

601 

328 

549 

324 

t .94 

538 

297 

535 

251 

503 

286 

392 

614 

358 

608 

313 

513 

244 

346 

572 

259 

543 

269 

498 

269 


583 368 

549 301 

251 -133 
671 399 

650 368 

560 137 

653 452 

669 381 

671 374 

593 141 

186 -18 
420 116 

221 -276 
652 378 

557 334 

589 377 

701 427 

564 -181 
678 329 

149 -86 

277 -259 
-19 -301 
693 475 

701 449 

762 485 

770 251 

707 469 

645 358 

675 362 

754 503 

78 . -251 
791 492 

840 544 

592 162 

789 256 

642 365 

749 542 

512 -227 
478 238 

441 -7 

163 -102 
733 420 

453 235 

222 -66 
571 213 

601 286 
514 -286 
94 -227 
776 -258 
233 -136 
518 47 

412 188 

689 329 

130 -237 
530 294 

594 374 

645 381 

656 427 

844 365 

570 309 


591 

360 

573 

313 

419 

63 

662 

325 

673 

-189 

739 

366 

709 

158 

706 

401 

763 

514 

511 

135 

206 

-8 

383 

-8 

233 

-194 

632 

342 

582 

255 

626 

426 

682 

286 

138 

-219 

671 

433 

474 

73 

292 

-145 

205 

-330 

722 

441 

685 

431 

696 

374 

484 

45 

789 

400 

625 

373 

667 

447 

709 

427 

526 

313 

717 

403 

792 

318 

530 

229 

475 

-387 

673 

336 

868 

545 

93 

-189 

501 

104 

511 

46 

632 

11 

803 

-495 

459 

117 

280 

-204 

695 

-448 

576 

221 

94 

-133 

50 

-156 

192 

-125 

251 

18 

589 

63 

999 

411 

712 

421 

-33 • 

-254 

575 

240 

598 

308 

639 

393 

635 

394 

607 

406 

513 

151 


601 

373 

632 

517 

317 

565 

535 

175 

490 

596 

331 

583 

242 

-17 

610 

725 

362 

717 

404 

-280 

253 

652 

382 

725 

757 

265 

519 

693 

412 

657 

739 

110 

453 

299 

66 

292 

226 

-11 

560 

637 

276 

557 

544 

240 

636 

648 

250 

647 

601 

175 

559 

251 

7 

443 

696 

283 

554 

363 

-419 

-17 

276 

-1 

227 

632 

231 

433 

742 

437 

730 

632 

425 

646 

591 

356 

639 

282 

-251 

255 

669 

334 

546 

652 

409 

634 

755 

349 

683 

703 

274 

516 

725 

314 

582 

727 

369 

603 

522 

-318 

110 

577 

217 

613 

232 

-347 

133 

889 

365 

669 

871 

415 

690 

518 

293 

594 

786 

299 

832 

895 

405 

877 

811 

390 

855 

141 

-296 

100 

691 

428 

764 

561 

197 

561 

103 

-312 

129 

721 

315 

636 

355 

36 

458 

329 

19 

585 

464 

211 

730 

528 

33 

449 

449 

205 

645 

942 

486 

757 

690 

378 

870 

551 

274 

549 

478 

234 

519 

581 

330 

739 

645 

377 

592 

640 

378 

653 

620 

387 

609 

506 

221 

582 
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347 

553 

-277 

162 

716 

34 

43 

511 

280 

-325 

-20 

-242 

367 

658 

408 

309 

0 

512 

26 0 


50 -204 
325 -249 
682 292 

508 254 

623 264 

571 346 


109 -256 
507 -15 

503 218 

478 192 

587 371 

725 261 


164 -65 

322 66 

705 197 

631 408 

589 285 

487 192 


313 49 

360 -158 
486 221 

716 350 

632 368 

593 293 


342 124 

83 -281 
458 202 

567 258 

599 318 

539 189 


365 

314 

608 

592 

589 

554 
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Table B-6. Output for Example 3 


MIN RANGE 1 200.000 
MAX RANGE=1494.000 
MIN MEAN=-160 . 000 
MAX MEAN® 753.000 
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MKM (2nd page) 


RANGE 


rainfloh CYCLES 


200.0 

250.0 

300.0 

350.0 

400.0 

450.0 

500.0 

550.0 

600.0 

650.0 

700.0 

750.0 

800.0 

850.0 

900.0 

950.0 
1000.0 

1050.0 

1100.0 

1150.0 

1200.0 

1250.0 

1300.0 

1350.0 

1400.0 

1450.0 

1500.0 


- 160.0 

0 

2 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


- 110.0 

5 

2 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


- 60.0 

5 

5 

5 

3 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


- 10.0 

2 

8 

3 

3 

1 

1 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


40.0 

2 

9 

2 

2 

2 

5 

2 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


90 


.0 

5 

9 

1 

3 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


140.0 

0 

5 

3 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 


190.0 

0 

5 

3 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

3 

2 

0 

0 

0 

0 

0 

1 

0 

0 

0 


TOTAL 

CYCLES 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
0 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 
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RAINFLOW CYCLES 


RANGE /xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxMEANxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


200.0 

250.0 

300.0 

350.0 

400 . 0 

450.0 

500.0 

550.0 

600.0 

650.0 

700.0 

750.0 

800.0 

850.0 

900.0 

950.0 
1000.0 

1050.0 

1100.0 

1150.0 

1200.0 

1250 . 0 

1300.0 

1350.0 

1400.0 

1450.0 

1500.0 


240.0 

3 

3 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

3 

4 

5 
3 
2 
3 
1 
2 
0 
0 
0 
0 
1 


290.0 

5 

6 
4 
3 

3 
2 
2 
0 
0 
0 
0 
0 
0 
1 
1 
2 
1 

4 
3 
0 
1 
0 
0 
0 
0 
0 
0 


340.0 

11 

13 

10 

6 

7 

0 

2 

4 

2 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


390.0 

24 

32 

15 

15 

10 

8 

2 

2 

2 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


440.0 

30 

37 

24 

17 

21 

9 

3 

6 

3 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


490.0 

33 

41 

23 

11 

9 

2 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


540.0 

25 

39 

20 

14 

7 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


590.0 

14 

14 

11 

6 

4 

2 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


"TOTAL 

CYCLES 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
0 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 
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RAINFLOH CYCLES 


RANGE 


200.0 

250.0 

300.0 

350.0 

400 . 0 

450.0 

500.0 

550.0 

600.0 

650.0 

700.0 

750.0 

800.0 

850.0 

900.0 

950.0 

1000.0 

1050.0 

1100.0 

1150.0 

1200.0 

1250.0 

1300.0 

1350.0 

1400.0 

1450.0 

1500.0 


640.0 

6 

10 

5 

3 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


690.0 

2 

2 

3 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


740.0 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


790.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


TOTAL 

CYCLES 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
Q 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 
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EXAMPLE 4 


The same history as Example 3 is analyzed using Option 1. Tables 
B-7 and B-8 show the input and output of this program, respectively. 

The output includes the filtered history as a peak/valley sequence and 
also the rain-flow cycles for the original history. 

In this case, filtering at a range of 400 reduced the length of the 
history from 1709 peaks and valleys to 361. The listing of the filtered 
history alternates between peaks and valleys, or valleys and peaks, and 
starts and ends with the largest absolute value. Note that range/mean 
matrix printed is the same as for Example 3. 
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Table B - 7 . Input for Example 4 


l 

400.0 
1709 
50 . 50 . 


0 

513 

292 

562 

267 

585 

31 $ 

52 $ 

299 

513 

299 

569 

272 

535 

298 

523 

291 

680 

$$$ 

683 

396 

678 

37 $ 

701 

276 

607 

331 

5$2 

339 

572 

274 

500 

285 

565 

27 $ 

$91 

237 

$96 

-76 

1$2 

-276 

71 

-309 

19 

-249 

137 

-93 

658 

275 

53 $ 

2$0 

612 

360 

690 

315 

689 

306 

698 

417 

666 

392 

699 

37 $ 

701 

368 

659 

3$5 

650 

$03 

610 

393 

637 

400 

618 

285 

630 

385 

616 

37 $ 

703 

- 17 $ 

136 

-331 

33 $ 

98 

510 

151 

542 

8 $ 

6$2 

398 

651 

-197 

$$ 

-197 

125 

-121 

228 

-17 

631 

387 

657 

$08 

623 

218 

553 

291 

6$0 

$08 

693 

$53 

657 

356 

642 

400 

698 

330 

571 

337 

689 

395 

6$7 

377 

59 $ 

227 

631 

299 

682 

411 

663 

$ 0 $ 

629 

306 

631 

306 

61 $ 

$09 

6 $8 

$ 1 $ 

683 

374 

575 

324 

691 

379 

585 

368 

621 

371 

610 

373 

582 

-8 

2$0 

-227 

248 

-205 

285 

-89 

$70 

2$8 

$60 

136 

519 

317 

517 

278 

503 

179 

615 

324 

559 

219 

679 

$21 

683 

- 25 $ 

8 $ 

-222 

2$3 

-19 

$95 

27 0 

629 

291 

56 0 

351 

6 57 

2$3 

608 

373 

663 

269 

653 

156 

663 

443 

690 

249 

612 

229 

602 

339 

738 

231 

$98 

237 

61 $ 

-167 

39 

-261 

437 

192 

503 

103 

682 

201 

589 

226 

$80 

135 

561 

3$7 

5$8 

149 

494 

260 

661 

158 

$67 

- 18 $ 

161 

-105 

555 

2$8 

$$9 

132 

712 

276 

529 

301 

52 $ 

309 

56 $ 

195 

$67 

35 

378 

- 29 $ 

188 

-308 

2$5 

-211 

607 

41 

526 

205 

506 

218 

513 

2$9 

603 

261 

712 

506 

727 

194 

400 

-244 

173 

-276 

610 

$05 

615 

360 

6$2 

$03 

6$1 

353 

556 

340 

655 

433 

669 

$22 

662 

$57 

796 

$03 

671 

78 

5$9 

285 

616 

376 

690 

18 

731 

108 

578 

282 

629 

170 

7 $8 

181 

578 

3$2 

60 $ 

261 

585 

330 

69 $ 

3$9 

61 $ 

313 

521 

290 

533 

266 

565 

291 

6$5 

-228 

1 

-286 

157 

-87 

725 

183 

71 $ 

393 

621 

19 $ 

$26 

152 

561 

235 

551 

213 

691 

136 

$87 

16 $ 

5$9 

2 $$ 

598 

162 

507 

2$9 

716 

339 

733 

194 

566 

-191 

120 

-179 

108 

-119 

181 

- 7 $ 

722 

$22 

639 

303 

671 

369 

715 

350 

605 

153 

765 

$09 

762 

355 

712 

398 

718 

362 

679 

448 

718 

$$6 

698 

3$2 

605 

188 

$98 

276 

537 

283 

522 

222 

546 

345 

57 0 

191 

500 

280 

658 

269 

685 

3$1 

551 

3 $$ 

618 

303 

612 

285 

666 

228 

53 $ 

303 

5$6 

-110 

110 

-362 

88 

-199 

173 

-130 

612 

333 

680 

3$1 

603 

283 

619 

30 

$06 

206 

$70 

260 

717 

218 

487 

190 

573 

270 

$76 

265 

607 

310 

78 $ 

573 

822 

$96 

825 

508 

887 

491 

769 

$86 

716 

3$2 

59 $ 

365 

578 

25 $ 

593 

355 

561 

324 

573 

153 

696 

$$3 

651 

372 

653 

-256 

39 

-173 

620 

213 

560 

217 

559 

274 

$91 

211 

571 

266 

6$0 

169 

572 

223 

$96 

221 

573 

264 

758 

358 

7 $3 

501 

777 

-161 

57 

- 1$6 

210 

-13 

216 

-15 

373 

133 

371 

151 

62 $ 

313 

628 

$05 

698 

-202 

1$6 

-92 

$83 

1$1 

551 

229 

551 

207 

533 

271 

55 $ 

218 

532 

210 

630 

179 

$89 

266 

618 

408 

650 

328 

599 

-211 

158 

- 20 $ 

68 

-179 

256 

-108 

11 $ 

- 11 $ 

$76 

228 

494 

178 

658 

218 

612 

2$5 

652 

$16 

658 

355 

678 

$28 

722 

437 

747 

542 

776 

-201 

$05 

188 

7 07 

$08 

6$0 

$22 

6$8 

$15 

727 

393 

610 

387 

722 

$17 

720 

$76 

688 

$ 8 $ 

691 

360 

570 

-350 

528 

216 

679 

475 

689 

27 $ 

562 

286 

653 

233 

$53 

253 

573 

282 

582 

356 

603 

234 

871 

$ 9 $ 

771 

511 

773 

$81 

727 

$91 

836 

566 

782 

398 

659 

329 

685 

387 

629 

218 

60 $ 

329 

628 

- 2$8 

27 $ 

-206 

710 

463 

663 

365 

695 

$06 

803 

$37 

739 

$$9 

679 

$52 

695 

$ 5 $ 

703 

26 0 

602 

-330 

3 $ 

-350 

188 

-180 

309 

-7 

6$8 

$17 

6$7 

$39 

683 

438 

699 

451 

712 

382 

698 

296 

625 

385 

680 

$05 

718 

$57 

673 

372 

744 

414 

618 

82 

597 

191 

570 

211 

522 

212 

517 

211 

$95 

232 

439 

235 

616 

256 

512 

27 0 

575 

31 $ 

791 

369 

787 

502 

769 

468 

785 

444 

798 

501 

732 

522 

723 

$62 

680 

- 26 $ 

68 

- 15 $ 

188 

-14 

242 

13 

582 

291 

618 

227 

537 

285 

555 

22 $ 

5$5 

329 

588 

374 

647 

242 

561 

353 

69 $ 

$33 

753 

-237 

313 

- 1$7 

226 

-163 

705 

329 

636 

376 

680 

383 

6$1 

$ 0 $ 

712 

367 

732 

290 

610 

318 

629 

315 

588 

353 

577 

265 

612 

368 

723 

-183 

267 

-208 

111 

-153 

$73 
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271 

576 

365 

618 

395 

605 

-283 

235 

-172 

687 

242 

521 

299 

532 

283 

596 

306 

629 

392 

639 

326 

583 

368 

591 

360 

601 

373 

632 

409 

690 

272 

496 

175 

470 

251 

549 

301 

573 

313 

517 

317 

565 

322 

529 

265 

586 

353 

555 

-215 

251 

-133 

419 

63 

535 

175 

490 

164 

517 

126 

56 0 

165 

726 

408 

671 

399 

662 

325 

596 

331 

583 

378 

749 

309 

635 

417 

677 

373 

650 

368 

673 

-189 

242 

-17 

610 

197 

548 

181 

523 

296 

582 

113 

560 

137 

739 

366 

725 

362 

717 

341 

607 

317 

604 

351 

567 

141 

653 

452 

709 

158 

404 

-280 

253 

-183 

199 

-28 

425 

167 

620 

400 

669 

381 

706 

401 

652 

382 

725 

365 

812 

395 

695 

416 

758 

469 

671 

374 

763 

514 

757 

265 

519 

27 

637 

212 

566 

275 

526 

325 

593 

141 

511 

135 

693 

412 

657 

281 

533 

-204 

199 

-185 

136 

-86 

186 

-18 

206 

-8 

739 

110 

453 

253 

506 

275 

511 

108 

516 

169 

420 

116 

383 

-8 

299 

66 

292 

-259 

-22 

-242 

244 

1 

222 

-17 

221 

-276 

233 

“194 

226 

-11 

560 

358 

567 

356 

599 

361 

607 

394 

652 

378 

632 

342 

637 

276 

557 

299 

583 

190 

506 

235 

57 0 

299 

557 

334 

582 

255 

544 

240 

636 

270 

619 

228 

583 

371 

593 

373 

589 

377 

626 

426 

648 

250 

647 

293 

717 

417 

701 

320 

709 

385 

701 

427 

682 

286 

601 

175 

559 

285 

533 

310 

535 

319 

555 

326 

564 

-181 

138 

-219 

251 

7 

443 

102 

490 

263 

534 

45 

693 

398 

678 

329 

671 

433 

696 

283 

554 

189 

489 

-367 

207 

-276 

-34 

-263 

149 

-86 

474 

73 

363 

-419 

-17 

-240 

-18 

-219 

181 

-470 

-19 

-293 

277 

-259 

292 

-145 

276 

-1 

227 

1 

566 

117 

516 

148 

415 

-406 

-19 

-301 

205 

-330 

632 

231 

433 

178 

582 

331 

835 

490 

748 

355 

693 

475 

722 

441 

742 

437 

730 

419 

722 

416 

769 

442 

775 

501 

701 

449 

685 

431 

632 

425 

646 

415 

683 

341 

769 

485 

738 

421 

762 

485 

696 

374 

591 

356 

639 

356 

675 

426 

705 

449 

688 

479 

770 

251 

484 

45 

282 

-251 

255 

-320 

438 

222 

639 

251 

647 

431 

707 

469 

789 

400 

669 

334 

546 

342 

705 

374 

684 

427 

698 

368 

645 

358 

625 

373 

652 

409 

634 

322 

624 

368 

659 

392 

748 

355 

675 

362 

667 

447 

755 

349 

683 

315 

693 

330 

607 

226 

744 

443 

754 

503 

709 

427 

703 

274 

516 

258 

588 

141 

457 

221 

428 

-164 

78 

-251 

526 

313 

725 

314 

582 

212 

733 

318 

629 

357 

752 

458 

791 

492 

717 

403 

727 

369 

603 

277 

636 

358 

562 

229 

770 

535 

840 

544 

792 

318 

522 

-318 

110 

-233 

288 

-100 

265 

-120 

703 

14 

592 

162 

530 

229 

577 

217 

613 

164 

593 

319 

632 

398 

623 

250 

789 

256 

475 

-387 

232 

-347 

133 

-141 

140 

-129 

544 

303 

644 

358 

642 

365 

673 

336 

889 

365 

669 

447 

813 

506 

817 

537 

776 

524 

749 

542 

868 

545 

871 

415 

690 

471 

789 

487 

706 

452 

669 

192 

512 

-227 

93 

-189 

518 

293 

594 

89 

576 

362 

789 

191 

449 

151 

478 

238 

501 

104 

786 

299 

832 

610 

829 

530 

734 

-318 

728 

83 

441 

-7 

511 

46 

895 

405 

877 

609 

876 

630 

888 

527 

797 

-239 

163 

-102 

632 

11 

811 

390 

855 

271 

578 

324 

675 

221 

632 

362 

733 

420 

803 

-495 

141 

-296 

100 

-140 

510 

251 

661 

366 

695 

161 

453 

235 

459 

117 

691 

428 

764 

-247 

217 

-301 

149 

-207 

157 

-74 

222 

-66 

280 

-204 

561 

197 

561 

153 

492 

27 0 

607 

240 

645 

202 

571 

213 

695 

-448 

103 

-312 

129 

-138 

508 

304 

603 

67 

626 

342 

601 

286 

576 

221 

721 

315 

636 

304 

669 

337 

605 

367 

567 

280 

514 

-286 

94 

-133 

355 

36 

458 

146 

506 

294 

577 

-11 

222 

-341 

94 

-227 

50 

-156 

329 

19 

555 

165 

398 

176 

822 

561 

801 

553 

776 

-258 

192 

-125 

464 

211 

730 

152 

834 

352 

642 

285 

587 

-222 

233 

-136 

251 

18 

528 

33 

449 

234 

626 

226 

480 

35 

594 

189 

518 

47 

589 

63 

449 

205 

645 

-229 

336 

87 

798 

529 

753 

51 

412 

188 

999 

411 

942 

486 

757 

394 

830 

464 

732 

330 

800 

452 

689 

329 

712 

421 

690 

378 

870 

462 

782 

116 

653 

-223 

54 

-223 

130 

-237 

-33 

-254 

551 

274 

549 

251 

658 

192 

640 

323 

578 

281 

530 

294 

575 

240 

478 

234 

519 

296 

645 

353 

664 

318 

575 

178 

594 

374 

598 

308 

581 

330 

739 

412 

615 

344 

601 

328 

549 

324 

645 

381 

639 

393 

645 

377 

592 

294 

538 

297 

535 

251 

503 

286 

656 

427 

635 

394 

640 

378 

653 

392 

614 

358 

608 

313 

513 

244 

844 

365 

607 

406 

620 

387 

609 
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346 

347 
162 

43 

-325 

367 

309 

0 


572 259 
553 -277 
716 34 

511 280 
-20 -242 
658 408 

512 260 


543 269 

50 -204 
325 -249 
682 292 

508 254 

623 264 

571 346 


498 269 

109 -256 
507 -15 

503 218 

478 192 

587 371 

725 261 


570 309 

164 -65 

322 66 

705 197 

631 408 

589 285 

487 192 


513 151 

313 49 

360 -158 
486 221 

716 350 

632 368 

593 293 


506 221 

342 124 

83 -281 
458 202 

567 258 

599 318 

539 189 


582 

365 

314 

608 

592 

589 

554 
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Table B-8. Output for Example 4 


MIN RANGE 5 200.000 
MAX RANGE=149A . 000 
MIN MEAN=-160 . 000 
MAX MEAN 5 753.000 
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FILTER HISTIRY-PEAK/VALLEY SEQUENCE 


FILTER LEVEL=400 . 000 


NUMBER OF POINTS IN FILTER HISTORY 3 361 
999.0 911.0 992.0 399.0 830.0 330.0 800.0 


870.0 

116.0 

653.0 

-254. 

664.0 

178.0 

739.0 

251. 

582.0 

-277.0 

716.0 

-249. 

705.0 

197.0 

608.0 

-325. 

593.0 

0.0 

701.0 

-309. 

703.0 

-331.0 

542.0 

84. 

698.0 

227.0 

691 .0 

-227. 

657.0 

243.0 

663.0 

156. 

682.0 

135.0 

561.0 

149. 

712.0 

-294.0 

188.0 

-308. 

727.0 

-244.0 

173.0 

-276. 

731.0 

108.0 

629.0 

170. 

725.0 

183.0 

714.0 

152. 

733.0 

-191.0 

722.0 

303. 

685.0 

-362.0 

680.0 

30. 

696.0 

-256.0 

620.0 

211. 

698.0 

-202.0 

630.0 

179. 

776.0 

-201.0 

727.0 

-350. 

628.0 

-248.0 

274.0 

-206. 

744.0 

82.0 

597.0 

191. 

647.0 

242.0 

753.0 

-237. 

723.0 

-183.0 

267.0 

-208. 

687.0 

175.0 

586.0 

-215. 

744.0 

-189.0 

610.0 

181 . 

739.0 

141.0 

709.0 

-280. 

637.0 

135.0 

693.0 

-204. 

516.0 

-259.0 

244.0 

-276. 

636.0 

228.0 

717.0 

-219. 

207.0 

-276.0 

474.0 

-419. 

292.0 

-145.0 

566.0 

-406. 

835.0 

355.0 

775.0 

341. 

255.0 

-320.0 

789.0 

322. 

725.0 

212.0 

733.0 

318. 

288.0 

-120.0 

703.0 

14. 

789.0 

-387.0 

232.0 

-347 . 

594.0 

89.0 

789.0 

104. 

728.0 

-7.0 

511.0 

46. 

632.0 

11.0 

811.0 

390. 

803.0 

-495.0 

141.0 

-296. 

217.0 

-301.0 

280.0 

-204, 

695.0 

-448.0 

103.0 

-312. 

721.0 

-286.0 

577 . 0 

-341. 

730.0 

152.0 

834.0 

-222. 

594.0 

999.0 

47.0 

589.0 

63. 


0 

65S.0 

192.0 

640.0 

0 

656.0 

244.0 

844.0 

0 

507.0 

-281.0 

682.0 

0 

716.0 

258.0 

725.0 

0 

658.0 

240.0 

701.0 

0 

651.0 

-197.0 

657.0 

0 

248 . 0 

-205.0 

683.0 

0 

690.0 

229.0 

738.0 

0 

661.0 

-184.0 

555.0 

0 

245.0 

-211.0 

607.0 

0 

796.0 

78.0 

690.0 

0 

748.0 

181.0 

694.0 

0 

691.0 

136.0 

598.0 

0 

715.0 

153.0 

765.0 

0 

717.0 

190.0 

887.0 

0 

640.0 

169.0 

777.0 

0 

650.0 

-211.0 

658.0 

0 

689.0 

233.0 

871.0 

0 

803.0 

-350.0 

712.0 

0 

791.0 

369.0 

798.0 

0 

313.0 

-163.0 

732.0 

0 

618.0 

-283.0 

235.0 

0 

535.0 

126.0 

726.0 

0 

582.0 

113.0 

560.0 

0 

253.0 

-183.0 

812.0 

0 

739.0 

110.0 

511.0 

0 

233.0 

-194.0 

652.0 

0 

534.0 

45.0 

696.0 

0 

181.0 

-470.0 

277.0 

0 

205.0 

-330.0 

632.0 

0 

769.0 

356.0 

770.0 

0 

755.0 

226.0 

754.0 

0 

791.0 

229.0 

840.0 

0 

592.0 

162.0 

613.0 

0 

889.0 

365.0 

871.0 

0 

786.0 

299.0 

832.0 

0 

895.0 

405.0 

888.0 

0 

855.0 

271.0 

675.0 

0 

695.0 

117.0 

764.0 

0 

561.0 

153.0 

645.0 

0 

603.0 

67.0 

626.0 

0 

585.0 

165.0 

822.0 

0 

528.0 

33.0 

626.0 

0 

645.0 

-229.0 

798.0 


329.0 

239.0 

151.0 

218.0 

192.0 

285.0 

218.0 
-259.0 
-261.0 

132.0 

91.0 

18.0 
-286.0 

162.0 

188.0 

153.0 
-161.0 

218.0 

218.0 

296.0 
-269.0 

265.0 
-172.0 

325.0 

137.0 

27.0 

108.0 

190.0 
-367.0 
-259.0 

178.0 
-251.0 
-251.0 
-318.0 

169.0 
-227.0 
-318.0 
-239.0 

221.0 

-297.0 

202.0 

221.0 

-258.0 

35.0 

51.0 


B— 35 



Table B-8 (3rd page) 


RAINFLOH CYCLES FOR ORIGINAL HISTORY 


NO RANGE LESS THAN OR EQUAL TO FILTER LEVEL=400 . GOO OCCUR IN FILTER HISTORY 

TOTAL 

RANGE /XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXMEANXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX CYCLES 


200.0 

250.0 

300.0 

350.0 

400.0 

450.0 

500.0 

550.0 

600.0 

650.0 

700.0 

750.0 

800.0 

850.0 

900.0 

950.0 

1000.0 

1050.0 

1100.0 

1150.0 

1200.0 

1250.0 

1300.0 

1350.0 

1400.0 

1450.0 

1500.0 


-160.0 

0 

2 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


- 110.0 
5 
2 
0 
1 
1 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


-6 0.0 
5 
5 
5 
3 
1 
1 
1 
1 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


- 10.0 

2 

8 

3 

3 

1 

1 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


40.0 

2 

9 

2 

2 

2 

5 

2 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


90.0 

5 

9 

1 

3 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


140.0 

0 

5 

3 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 


190.0 

0 

5 

3 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

3 

2 

0 

0 

0 

0 

0 

1 

0 

0 

0 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
0 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 


c-z- 
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NO 

RANGE 


RAINFLOW CYCLES FOR ORIGINAL HISTORY 
RANCH LESS THAN OR EQUAL TO FILTER LEVEl=AOO . OOO OCCUR IN FILTER 


200.0 

250.0 

300.0 

350.0 

400.0 

450.0 

500.0 

550.0 

600.0 

650.0 

700.0 

750.0 

800.0 

850.0 

900.0 

950.0 

1000.0 

1050.0 

1100.0 

1150.0 

1200.0 

1250.0 

1300.0 

1350.0 

1400.0 

1450.0 

1500.0 


240.0 

3 

3 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

3 

4 

5 
3 
2 
3 
1 
2 
0 
0 
0 
0 
1 


290.0 

5 

6 
4 
3 
3 
2 
2 
0 
0 
0 
0 
0 
0 
1 
1 
2 
1 
4 
3 
0 
1 
0 
0 

0 

0 

0 

0 


340.0 

11 

13 

10 

6 

7 

0 

2 

4 

2 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


390.0 

24 

32 

15 

15 

10 

8 

2 

2 

2 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


440.0 

490.0 

540.0 

30 

33 

25 

37 

41 

39 

24 

23 

20 

17 

11 

14 

21 

9 

7 

9 

2 

2 

3 

2 

1 


6 

3 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


*xxxxxxxxx 


590.0 

14 

14 

11 

6 

4 

2 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


HISTORY 

total 

CYCLES 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
0 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 
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Table B-8 (5th page) 


RAINFLOW CYCLES FOR ORIGINAL HISTORY 


NO RANGE LESS THAN OR EQUAL TO FILTER LEVEL=400 . OOO OCCUR IN FILTER HISTORY 

RANGE J?JlES 


640.0 690.0 740.0 790.0 

200.0 ’ 


250.0 10 

300.0 5 

350.0 3 

400.0 0 

450.0 0 

500.0 2 

550 . 0 0 

600.0 0 

650.0 0 

700.0 0 

750.0 0 

800.0 0 

850.0 0 

900.0 0 

950.0 0 

1000.0 0 

1050.0 0 

1100.0 0 

1150.0 0 

1200.0 0 

1250.0 0 

1300.0 0 

1350.0 0 

1400.0 0 

1450.0 0 

1500.0 0 


2 2 0 

3 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

10 0 
0 0 0 

0 0 0 

ooo 

0 0 0 

ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
ooo 
0 0 0 
ooo 
ooo 
ooo 


173 

244 

134 

91 

68 

33 

21 

18 

10 

3 

3 

4 
0 
8 

4 
7 
9 
9 

5 
3 
3 
2 
0 
1 
0 
0 
1 
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EXAMPLE 5 


A history containing 1021 peak/valley points is used. This is the 
Filtered version (filter level = .46) of a simulated helicopter combat 
maneuver loading history. Option 4 is used, so that the results are 
given in the form of a compact 32 x 32 matrix containing the maximum and 
minimum (peak and valley) values of the rain-flow cycles. Tables B-9 
and B-10 show the input and output of this program, respectively. Note 
that the history is converted using linear interpolation to have a 
minimum value of 1 and a maximum value of 32. 
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Table B-9 


. Input for Example 5 


4 

1021 

1.000 0.430 

0.919 0.467 

0.802 0.135 

0.835 0.380 

0.836 - 0.258 
0.287 - 0.216 
0.420 - 0.092 
0.406 - 0.084 
0.804 0.335 

0.820 0.264 

0.765 0.200 

0.894 0.423 

0.932 0.395 

0.514 0.018 

0.545 0.073 

0 . 866 0.304 

0.834 0.363 

0.837 0.334 

0.727 0.192 

0.706 0.191 

0.616 0.116 
0.598 0.080 

0.860 - 0.035 
0.696 0.157 

0.668 - 0.071 
0.436 - 0.146 
0.721 0.254 

0.799 0.303 

0.742 0.167 

0.725 0.116 

0.731 0.051 

0.739 0.079 

0.705 0.109 

0.708 0.220 

0.939 0.398 

0.995 0.365 

0.872 0.362 

0.826 0.371 

0.811 0.039 

0.611 0.010 
0.287 - 0.201 
0.364 - 0.185 
0.304 - 0.155 
0.301 - 0.266 
0.742 0.271 

0.703 0.190 

0.624 0.099 

0.577 0.084 

0.337 - 0.124 
0.540 - 0.066 
0.497 - 0.060 


0.973 

0.957 

0.925 

0.992 

0.721 

0.159 

0.896 

0.397 

0.933 

- 0.033 

0.355 

- 0.222 

0.388 

- 0.067 

0.393 

- 0,067 

0.891 

0.306 

0.758 

- 0.010 

0.697 

0.230 

0.983 

0.923 

0.936 

0.380 

0.511 

0.096 

0.726 

0.235 

0.791 

0.329 

0.857 

0.392 

0.876 

0.157 

0.669 

0.191 

0.712 

0.183 

0.621 

0.088 

0.550 

0,073 

0.510 

0.032 

0.719 

0.059 

0.630 

- 0 . 098 

0.903 

- 0 .129 

0.722 

0.219 

0.839 

0.296 

0.797 

0.183 

0.700 

0.189 

0.732 

0.129 

0.711 

0 . 063 

0.785 

0.075 

0.826 

0.361 

0.882 

0.365 

0.979 

0.382 

0.839 

0.381 

0.839 

0.392 

0.583 

0.059 

0.592 

0.000 

0.268 

- 0.219 

0.316 

- 0.135 

0.322 

- 0.153 

0.659 

0.182 

0.776 

0.259 

0.692 

0,197 

0.603 

0.197 

0,591 

0.015 

0.930 

- 0.081 

0.561 

- 0.191 

0.952 

- 0.089 


0,978 

0.956 

0.971 

0.926 

0.691 

- 0.091 

0.860 

0.369 

0.922 

- 0.033 

0.558 

0.059 

0.999 

- 0.139 

0.387 

- 0.119 

0.805 

0.390 

0.793 

0.289 

0.771 

0.321 

0.909 

0.369 

0.867 

0.398 

0.559 

- 0.516 

0.733 

0.260 

0.319 

0.398 

0.969 

0.350 

0.710 

0.199 

0.785 

0.163 

0.690 

0.165 

0.627 

0.165 

0.569 

0.091 

0.995 

- 0.201 

0.657 

0,132 

0.626 

- 0.091 

0.339 

- 0.171 

0.828 

0.335 

0.788 

0.309 

0.781 

0.259 

0.692 

0.076 

0.656 

0.118 

0.698 

0.153 

0.561 

0.098 

0.832 

0.373 

0.911 

0.992 

0.976 

0.371 

0.882 

0.919 

0.896 

0.337 

0.535 

- 0.097 

0.951 

- 0.299 

0.275 

- 0.190 

0.36 0 

- 0.103 

0.336 

- 0.137 

0.702 

0.220 

0.783 

0.279 

0.630 

0.106 

0.710 

0.092 

0.972 

0.021 

0.912 

-0 . 125 

0 . 387 

- 0.072 

0.911 

- 0.099 


0,939 

0.467 

0.889 

0.357 

0.799 

0.001 

0,832 

0.365 

0.939 

- 0.174 

0.528 

- 0.055 

0.381 

- 0.121 

0.371 

- 0.087 

0.796 

0.316 

0,797 

0.333 

0.886 

0.335 

0.860 

0.383 

0.806 

0.197 

0.381 

- 0.215 

0.776 

0.294 

0.800 

0.301 

0.909 

0.254 

0.776 

0.209 

0.775 

0.179 

0.679 

0.149 

0.635 

0.122 

0.625 

- 0.016 

0.603 

0.009 

0.727 

0.109 

0.527 

- 0.112 

0.280 

- 0.382 

0.898 

0.365 

0.780 

0.303 

0.760 

0.099 

0,619 

0.128 

0.755 

0.079 

0.619 

0.144 

0.691 

0.146 

0.832 

0.375 

0.997 

0.393 

0.978 

0.412 

0.889 

0.346 

0.790 

0.321 

0.555 

0.005 

0.277 

- 0.237 

0.280 

- 0.174 

0.352 

- 0.115 

0.327 

- 0.142 

0.733 

0.231 

0.757 

0.270 

0.629 

0.165 

0.711 

0.047 

0.512 

- 0.199 

0.917 

- 0.074 

0.382 

- 0.082 

0.930 

- 0.131 


0.929 

0.909 

0.887 

0.339 

0.750 

0.127 

0.899 

0.350 

0.293 

- 0.180 

0.398 

- 0.096 

0.356 

- 0.095 

0.379 

- 0.129 

0.835 

0.329 

0.803 

0.199 

0.879 

0.373 

0.879 

0.922 

0.676 

0.002 

0,376 

- 0.109 

0.819 

0.317 

0.828 

0.303 

0.863 

0.350 

0.709 

0.223 

0.701 

0.216 

0.631 

0.127 

0.581 

0.058 

0.596 

- 0.237 

0.631 

0.066 

0.589 

- 0.011 

0.529 

- 0.109 

0.298 

- 0.193 

0.871 

0.232 

0.775 

0.275 

0.692 

0.108 

0.600 

0.088 

0.693 

0.053 

0.681 

0.128 

0.771 

0.230 

0.889 

0.369 

0.901 

0.373 

0.883 

0.390 

0.881 

0.398 

0.885 

0.287 

0.958 

0.001 

0.303 

- 0.215 

0.312 

- 0.163 

0,352 

- 0.198 

0.332 

- 0.169 

0.796 

0.237 

0.729 

0.196 

0.730 

0.139 

0.566 

0,115 

0.321 

- 0.158 

0.992 

- 0.026 

0.968 

- 0.059 

0.382 

- 0.131 
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table B-9 (2nd page) 


0.356 

- 0.099 

0.405 

- 0.116 

0.488 

- 0.144 

0.460 

- 0.086 

0.440 

- 0.121 

0.553 

- 0.115 

0.501 

- 0.077 

0.487 

- 0.041 

0.495 

- 0.012 

0.470 

- 0.026 

0.519 

- 0.019 

0.532 

- 0.013 

0.459 

- 0.022 

0.542 

- 0.062 

0.575 

- 0.016 

0.549 

- 0.017 

0.542 

0.024 

0.573 

0.012 

0.681 

0.053 

0.640 

0.108 

0.595 

0.139 

0.602 

0.028 

0.613 

0.015 

0.635 

- 0.079 

0.582 

0.013 

0 . 530 

0.052 

0.620 

0.161 

0.647 

0.180 

0.635 

0.164 

0.685 

0.180 

0.709 

0.129 

0.701 

0.127 

0.643 

0.104 

0.615 

- 0.075 

0.471 

- 0 . 049 

0.432 

- 0 . 032 

0.423 

- 0.060 

0.396 

- 0.095 

0.458 

- 0.055 

0.429 

- 0.083 

0 . 494 

- 0.080 

0.486 

- 0.040 

0.515 

- 0.047 

0.505 

- 0.101 

0.477 

- 0.049 

0 . 509 

- 0 . 065 

0.541 

- 0.087 

0.590 

- 0.084 

0.579 

0.077 

0.595 

- 0.098 

0 .654 

- 0.133 

0.612 

- 0.087 

0.618 

- 0.137 

0.572 

- 0.134 

0.490 

- 0.043 

0.435 

- 0.107 

0.366 

- 0.099 

0.536 

0.000 

0.456 

0.005 

0.695 

0.042 

0.648 

0.117 

0.731 

0.238 

0.803 

0.275 

0.811 

0.193 

0.825 

0.248 

0.700 

0.171 

0.887 

0.217 

0.762 

0.185 

0.721 

0.236 

0.701 

0.200 

0.732 

0.227 

0.873 

0.197 

0.800 

0.192 

0.684 

0.221 

0.738 

0.194 

0.742 

0.161 

0.750 

0.223 

0.718 

0.181 

0.742 

0.253 

0.732 

0.241 

0.731 

0.194 

0.730 

0.228 

0.776 

0.265 

0.790 

0.226 

0.795 

0.238 

0.717 

0.244 

0.765 

0.270 

0.768 

0.226 

0.732 

0.256 

0.748 

0.216 

0.755 

0.286 

0.787 

0.272 

0.727 

0.232 

0.738 

0.187 

0.749 

- 0.341 

0.299 

- 0.155 

0.311 

- 0.180 

0.650 

- 0.198 

0.318 

- 0.182 

0.281 

- 0.231 

0 . 442 

- 0 . 073 

0.453 

- 0.085 

0.471 

- 0.061 

0.467 

0.000 

0.483 

0.002 

0.531 

0.044 

0.546 

- 0.041 

0.541 

0.010 

0.564 

0.047 

0.594 

0.040 

0.534 

0.033 

0.570 

0.059 

0.556 

0.101 

0.562 

0.097 

0.587 

0.110 

0.601 

0.012 

0.625 

- 0.259 

0.683 

0.220 

0.748 

0.287 

0.849 

0.328 

0.858 

0.401 

0.870 

0.389 

0.892 

0.375 

0.900 

0.446 

0.917 

0.436 

0.948 

0.475 

0.971 

0.403 

0.951 

0.460 

0.951 

0.404 

0.884 

0.365 

0.975 

0.447 

0.950 

0.440 

0.914 

0.449 

0.922 

0.414 

0.926 

0.456 

0.941 

0.450 

0.910 

0.408 

0.880 

0.352 

0.824 

0.342 

0.865 

0.358 

0.864 

0.381 

0.936 

0.346 

0.872 

- 0.213 

0.559 

0.061 

0.513 

0.040 

0 . 509 

- 0.020 

0.466 

- 0.151 

0.775 

0.305 

0.821 

0.229 

0.788 

0.333 

0.834 

0.333 

0.819 

0.362 

0.829 

0.302 

0.833 

0.324 

0.846 

0.290 

0.802 

0.252 

0.813 

0.309 

0.775 

0.175 

0.797 

0.239 

0.780 

0.273 

0.737 

0.276 

0.789 

0.304 

0.846 

0.291 

0.742 

0.285 

0.767 

0.240 

0.775 

0.211 

0.722 

0.221 

0.698 

0.205 

0.707 

0.174 

0.690 

0.002 

0.673 

0.172 

0.790 

0.225 

0.760 

0.294 

0.795 

0.308 

0.789 

0.334 

0.822 

0.351 

0.811 

0.303 

0.802 

0.289 

0.782 

0.301 

0.784 

0.325 

0.783 

0.299 

0.847 

0.324 

0.822 

0.343 

0.896 

0.373 

0.875 

0.394 

0.943 

0.383 

0.857 

0.350 

0.924 

0.372 

0.894 

0.421 

0.874 

0.383 

0.878 

0.357 

0.846 

0.379 

0.908 

0.377 

0.913 

0.434 

0.887 

0.390 

0.879 

0.379 

0.878 

0.388 

0.839 

0.334 

0.842 

0.381 

0.896 

0.395 

0.901 

0.393 

0.872 

0.374 

0.946 

0.398 

0.881 

0.392 

0.897 

0.338 

0.958 

0.402 

0.912 

0.337 

0.877 

0.398 

0.878 

0.379 

0.838 

0.380 

0.920 

0.389 

0.856 

0.383 

0.858 

0.377 

0.910 

0.398 

0.860 

0.399 

0.881 

0.428 

0.901 

0.362 

0.910 

0.417 

0.928 

0.414 

0.897 

0.394 

0.956 

0.422 

0.943 

0.407 

0.949 

0.425 

0.954 

0.394 

0.922 

0.423 

0.955 

0.337 

0.905 

0.444 

0.934 

0.413 

0.932 

0.405 

0.926 

0.383 

0.907 

0.362 

0.860 

0.004 

0.494 

- 0.016 

0.771 

0,141 

0.900 

0.431 

0.952 

0.406 

0.961 

0.470 

0.952 

0.440 

0.995 

0.401 

0.944 

0.465 

0.972 

0.414 

0.952 

0.485 

0.958 

0.437 

0.987 

0.399 

0.990 

0.414 

0.980 

0.457 

0.981 

0.407 

0.962 

0.452 

0.959 

0.407 

0.910 

0.437 

0.966 

0.440 

0.945 

0.453 

0.932 

0.421 

0.939 

0.438 

0.920 

0.000 


1 .000 
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Table B-10. Output for Example 5 


RAINFLOW CYCLES 


PEAK 


1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8 . 0 
9.0 

10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

17.0 

18.0 

19.0 

20.0 
21.0 
22.0 

23.0 

24.0 

25.0 

26.0 

27.0 

23.0 

29.0 

30.0 

31.0 

32.0 


1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


2.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


3.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


4.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 



TOTAL 

CYCLES 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

17 

18 
20 
22 
27 
33 
33 
24 
48 
47 
51 
50 
60 
41 
11 
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5 Tabte B-10 (2nd page) 


PEAK 


1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

17.0 

18.0 

19.0 

20.0 
21.0 
22.0 

23.0 

24.0 

25.0 

26.0 

27.0 

28.0 

29.0 

30.0 

31.0 

32.0 


9.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

6 

5 

1 

4 

4 

1 

1 

1 

0 

0 

0 

0 

0 

0 


RAIN FLOW CYCLES 


10.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

9 

9 

5 

3 

5 

1 

0 

0 

0 

0 

0 

0 

0 


11.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

6 
6 
3 
2 
1 
0 
0 
0 
1 
1 
1 
0 


12.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

8 

9 

4 

1 

1 

1 

0 

1 

0 

0 

1 


13.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

8 

4 

3 

4 
2 
1 
0 
1 
0 
0 


14.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 
9 

6 
8 
3 
1 
0 
0 
0 
0 


TdTAL 

CYCLES 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

17 

18 
20 
22 
27 
33 
33 
24 
48 
47 
51 
50 
60 
41 
11 


B-43 



Table B-10 (3rd page) 


RAINFLOH CYCLES 


PEAK 


1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

17.0 

18.0 

19.0 

20.0 
21.0 
22.0 

23.0 

24.0 

25.0 

26.0 

27.0 

28.0 

29.0 

30.0 

31.0 

32.0 


17.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

16 

6 
1 
0 
0 
0 


18.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

22 

8 
2 
4 
0 


19.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

25 

19 

6 

2 


20.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

25 

16 

4 


21.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

12 

1 


22.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


23.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


24.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


TOTAL 

CYCLES 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

17 

18 
20 
22 
27 
33 
33 
24 
48 
47 
51 
50 
60 
41 
11 


B-44 



Table B-10 (4th page) 

“ T/y-' .• > 

RAINFLOH CYCLES 

PEAK /xxxxxxxxxxxxxxxxxxxxxxxxxxxxxvAL L EYxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


TOTAL 

CYCLES 


1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 
9.0 

10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

17.0 

18.0 

19.0 

20.0 
21.0 
22.0 

23.0 

24.0 

25.0 

26.0 

27.0 

28.0 

29.0 

30.0 

31.0 

32.0 


25.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


26.0 27.0 28.0 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0-0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


29.0 30.0 31.0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 _ 0 

0 0 0 

0 0 0 

0 0 0 

o 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 '0 0 

0 0 0 

0 0 0 


32.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

17 

18 
20 
22 
27 
33 
33 
24 
48 
47 
51 
50 
60 
41 
11 


B-45 



NASA 

Nol*>V3< fi^rcT^ifes and 
‘^rxe Aonwny jo* 


Report Documentation Page 


1. Report No. 


NASA CR-181941 


2. Government Accession No. 


4. Title and Subtitle 

Fatigue Life Estimates for Helicopter Loading Spectra 


7. Authorfs) 

A.K. Khosrovaneh, N.E. Dowling, 
A.P. Berens, and J.P. Gallagher 


9. Performing Organization Name and Address 

Virginia Polytechnic Institute and State University 
Engineering Science and Mechanics Department 
Blacksburg, Virginia 24061 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


3, Recipient's Catalog No. 


5. Report Date 


December 19K9 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

505-63-01-05 


11. Contract or Grant No. 

NAG1-&22 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Langley Technical Monitor: R. A. Everett, Jr. 

A.K. Khosrovaneh and N.E. Dowling: Virginia Polytechnic Institute and State University, 
Blacksburg, V A 

A.P. Berens and J.P. Gallagher: University of Dayton Research Institute, Dayton, OH 


16. Abstract 

Helicopter loading histories applied to notched metal samples are used as examples, and their fatigue lives 
are calculated by using a simplified version of the local strain approach. This simplified method has the 
advantage that it requires knowing the loading history in only the reduced form of ranges and means and 
number of cycles from the rain-flow cycle counting method. The calculated lives compare favorably with 
test data. 


17. Key Words (Suggested by Authorlsi) 

Rain-flow 
Helicopter 
Local Strain 
Fatigue 


18. Distribution Statement 


Unclassified - Unlimited 
Subject Category - 39 


V^UillUIUllVV UUIMU^V 

19. Security Classif. (of this report) 

20. Security Classif (of this page) 

21. No of pages 

22. Price 

Unclassified 

Unclassified 

107 

AQ6 


NASA FORM 1«26 OCT 86 




